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EDITORIAL 


A Técnica conseguiu finalmente recobrar tempo perdido. 
Este terceiro número, ainda que ligeiramente atrasado, ainda vai a tempo. 


Tece-se o quarto número da Técnica pela singularidade de ser um número 
especial, apontando o futuro. Isso tem-nos obrigado a esforços grandes mas, 


esperamos, os quatro números de 1991 sairão em 1991. 
Problemas materiais alguns. 


As melhores relações com o Conselho Directivo do Instituto Superior 


Técnico, é grato dizê-lo. 


E como mais de uma vez o temos feito, um convite à generalizada 


comparticipação de todos os docentes do IST. 
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Normas de Admissão de Trabalhos e 
Instruções para os Autores 


A revista TECNICA, em princípio, publicará 
anualmente quatro números de natureza Lécnico- 
cientifica. 


Atendendo do contexto onde se insere e público a 
que, em conformidade e de modo principal, se 
destina, a revista integra-se a partir das três 
seguintes formas de publicação: 


ARTIGOS POR CONVITE 

São da iniciativa da Direcção da Técnica é 
destinam-se, em princípio, a trabalhos de sintese no 
domínio das ciências básicas ou ciências da 
engenharia Além do mais, procura-se que estes 
artigos se constituam fontes recomendadas de 
estudo aos ulunos e de consulta para 
investigadores « docentes. 


ARTIGOS TECNICO-CIENTIFICOS 

São artigos submetidos a publicação. Destinam-se 
essencialmente a dar noticia sobre aquisições 
actuais da ciência e da técnica mas em Lermos que 
se vá além da simples divulgação, sendo desejável 
que nesses trabalhos figurem parâmetros e 
referências que permitam elaborar juízos críticos € 
comparativos. Hundamentalmente, uma noticia 
sucinta, com carácter crítico, dando conta de uma 
aquisição recente e importante para a técnica é a 
marca [fundamental destes artigos que vão 
traduzir, em princípio, investigação científica 
realizada pelo autor 


NOTAS CIENTIFICAS 
Destinam-se essencialmente à apresentação de 
novos resultados no domínio da investigação 
científica ou do desenvolvimento tecnológico. 
Pretende-se uma nota curta, realçando os 
resultados alcançados e referenciando 
convenientemente trabalhos afins, evitando a 
exposição longa de antecedentes, que deverá ser 
remetida pura referências. Deste modo, a 
originalidade que se pretende e se traduzirá pelos 
resultados novos que o autor apresenta é a marca 
fundamental de uma nota científica. 


Sempre que tal se ofereça conveniente, a TÉCNICA 
ouvirá especialistas sobre a essência e estrutura dos 
trabalhos submetidos a publicação, com natural 
excepção dos Artigos por Convite. Quando o autor ou 
autores forem alunos é necessário que do trabalho 
conste a indicação de um supervisor de reconhecida 
idoneidade cientifica e técnica que tenha 
comparticipado na análise crítica e prévia discussão 
do trabalho em causa 


Os Artigos por Convite e os Artigos Tecenico- 
Científicos, deverão ser escritos em língua 
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portuguesa. As Notas Científicas poderão ser 
escritas em inglês ou francês. 


Os autores deverão atender na elaboração dos 
arligos às seguintes condições: 

-os Artigos por Convite não poderão exceder 30 
páginas dactilografadas, incluindo o espaço 
reservado a Ilustrações; 

-os Artigos Tecnico-Cientificos não poderão exceder 
15 páginas dactilografadas, incluindo o espaço 
reservado a ilustrações; 

-as Notas Científicas não poderão exceder 8 páginas 
dactilografadas, incluindo o espaço reservado a 
ilustrações; 

-nome exacto com que desejam figurar; 

-ordem de entrada dos autores, havendo mais do que 
um; 

-instituição(ões) a que pertençam e que desejem 
referir, indicando claramente o nome dessa 
instituição, morada e sigla. 


Os trabalhos devem ser dactilografados a dois 
espaços, em folhas modelo A4, devidamente 
numeradas. 


Os trabalhos devem ser enviados em triplicado, 
recomendando-se aos autores que conservem uma 
cópia. 


Conjuntamente com o trabalho, devem ser 

enviados: 

—título do trabalho em português e inglês; 

-resumo em português e abstract em inglês, pelo 
menos; 

-nome, direcção e telefone do autor e demais dados 
que facilitem o processo de revisão de provas. 


Us autores deverão assegurar-se, antes de 
enviarem os originais, que estes não contenham 
erros, falta de dados, ete., pois que, na fase de revisão 
de provas, já não se admite inclusão de novo original 
e apenas se corrigirão erros de composição 
tipográfica, 


As referências bibliográficas serão indicadas no 
texto por ordem numérica, usando algarismos 
árabes entre colchetes [|], e a lista das mesmas, no 
final do artigo, obedecerá às seguintes normas; 


Artigos de revistas: 
Nome(s) do(s) autor(es) e iniciais (ano de 
publicação entre parêntesis). Nome do artigo, 
Nome da revista, número do volume, página(s), 


Livros: 

Nomels) do(s) autor(es) e iniciais (ano de 
publicação entre parêntesis). Título, editora, 
página(s). 
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Convem distinguir claramente entre referências e 
outra bibliografia que se recomende. 

Trabalhos em que a matéria apresentada não esteja 
referenciada não são, em princípio, publicados. Este 
é um ponto fundamental em que se insiste. 


Os quadros ou tabelas, numerados em romano, 
devem figurar no texto dactilografado, podendo ser 
encimados por um título conciso. 


Os gráficos devem, de preferência, ser desenhados a 
tinta da china, em folhas separadas e de modo a 
estarem prontos para serem reproduzidos com a 
maior nitidez e qualidade gráfica possível, 
Quaisquer letras ou algarismos que neles figurem 
deverão ser suficientemente grandes para permitir a 


sua redução que, em princípio, fica ao critério do 
editor. 

As folhas com os gráficos, figuras ou fotografias, 
deverão ser numeradas, indicando-se o trabalho a 
que dizem respeito. Além disso, no texto, indicar-se- 
ão local para a respectiva inclusão, 


NOTA: A TECNICA salienta a necessidade de os 
autores, quando for o caso, obterem em tempo as devidas 
autorizações de reprodução de excertos de outras obras, 
sendo esta matéria da inteira responsabilidade do autor. 


A TÉCNICA oferecerá cinco exemplares da Revista, 
por cada trabalho publicado, independentemente do 
número de autores. 
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RESUMO 


Vivemos um tempo de mudança. Assistimos ao 
"desmoronamento" de um mundo e ao construir dos 
alicerces dum novo. É a revolução Pós-Industrial, cuja 
implantação continuará a provocar grandes 
convulções e a obrigar a uma nova forma de olhar o 
mundo: global e intregada (visão sistémica). 

Esta mudança compota em si alterações e desafios a 
que os recursos humanos e organizacionais se terão 
que adaptar. As Universidades e as suas escolas, 
sobretudo as técnicas, não podem ficar indiferentes a 
esta revolução, nem ao que ela implica ao nível dos 
perfis curriculares. 

O IST, a maior Escola de Engenharia do País, tem 
responsabilidades acrescidas e, o Lipo de respostas que 
encontrar, para vencer os novos desafios, têm uma 
enorme relevância no contexto de funcionamento e 
necessidades das organizações industrias. 


] Mudam-se os tempos... 

As mudanças tecnológicas e sociais são uma 
constante, ao longo dos tempos e, a velocidade moderada 
dessas alterações, foi permitindo, às pessoas e às 
organizações, o ajustamento e a adaplação a essas 
mudanças. 

No entanto, desde o início do século XX, essa 
velocidade foi progressivamente aumentando, vindo a 
produzir-se uma variação muito significativa a partir 
dos finais dos anos 40, conduzindo 4 uma nova revolução 
- à Revolução Pós-Industrial que está em processo 
profundo de implementação. 

Esta revolução, pois trata-se de uma autêntica 
revolução, baseia-se na disponibilidade e utilização de 
três tecnologias fundamentais: 


- Desenvolvimento de instrumentos geradores de 
símbolos que representam as propriedades de 
objectos ou acontecimentos (data) 

- Transmissão de símbolos (tecnologias de 
comunicação) 

- Manipulação lógica de símbolos (computadores) 


Com a utilização simultânea destas três tecnologias, 
foi possivel utilizar sistemas para controlar outros 
sistemas, para automatizar. E aqui está a diferença 
fundamental entre a Revolução Industrial e a Revolução 
Pós-Industrial. Enquanto a primeira tinha a ver com a 
substituição do músculo - do homem pela máquina 


Original submetido para publicação em 22 de Outubro de 1991 


ABSTRACT 


We are living a time of changes. We are witnessing 
the "collapse” of a world and the building of another 
world's foundations. It is the Post-Industrial 
Revolution the implantation of which will go on 
causing upheavals and originate a new way of looking 
at the world: a global and integrated one (systems 
view). 

This change implies alterations and challenges to 
which the human and organizational resources will 
have to adapt. The Universities and their schools, 
mainly the technical ones, cannot remain indifferent 
to this revolution and to what it means in terms of 
curricula, 

The IST, the country's largest school of Engineering 
has additional responsabilities and the answers it 
provides to overcome Lhe new challenges will have an 
enormous importance in the performance and needs of 
the industrial organizations. 


como fonte de trabalho, fazendo tarefas repetitivas 
muito simples, produzindo a alienação do trabalho - a 
segunda tem a ver com a substituição da mente, 
permitindo a progressiva subsituição do homem nas 
tarefas rotineiras e repetitivas - automação. 

Esta revolução foi acompanhada, como sempre 
acontece, por uma forma diferente/própria de olhar e 
interpretar o mundo. 

Enquanto que durante a “Idade da Máquina” 
(anterior à Revolução Pós-Industrial) a compreensão 
(visão do mundo) era o resultado da compreensão das 
partes suas constituintes - doutrina determinista e 
método analítico, com a 2º Guerra Mundial, que 
"obriga" os cientistas a sair dos seus laboratórios e a 
embrenhar-se no “mundo real" para resolver problemas 
de grande dimensão em grupos interdisciplinares, e com 
base nalguns pressupostos que vinham a ser elaborados 
desde o início do presente século, começa a tomar forma 
uma visão sistémica do mundo - Idade dos Sistemas. 

Esta visão, cada vez mais prevalecente, do Mundo é 
caracterizada por fazer preceder a análise pela síntese, 
procurando primeiro identificar e explicar o supra 
sistema, de que o sistema a estudar faz parte, e, só 
depois, procurar explicar este. 

A visão sistémica faz assim intervir como elemento 
determinante no funcionamento dum qualquer sistema o 
ambiente em que este se insere, e só depois identifica os 
subsistemas constituintes e as suas interrelações. 

Em resumo, estamos em plena revolução, onde se olha 
para o mundo de uma forma global e integrada - visão 
sistémica, e se dispõe dum conjunto de tecnologias que 
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permitem a progressiva substituição do homem nas 
Larefas rotineiras € repetitivas. 

Esta nova revolução, com uma elevada Laxa de 
aceleração da mudança está a pôr em causa “velhas” e 
assumidas convicções, modos de estar, de produzir, etc... 

Assistimos ao “desmoronamento” dum mundo e ao 
construir dos alicerces dum novo. Este parto não está a 
ser feito sem dôr, e, como em qualquer revolução, 
continuará a provocar grandes convulsões, 

Gerações (humanas e organizacionais de todos os 
sectores) estão lutando desesperadamente para se 
adaptarem às novas circunstâncias garantindo a 
sobrevivência. Nem todas o conseguirão ! 


2 Mudam-se as organizações 

Às empresas, à luz du visão sistêmica, são 
organizações que integram em si três níveis de 
objectivos: 


-os do Sistema (suprasistema) em que estão 
englobados - (sociais/ambientais) 

- os do Sistema em análise (organizacionais) 

-o das pessoas que são parte desse sistema 
(individuais) 


Por outro lado as mudanças aceleradas exigem das 
empresas uma atenção constante ao ambiente em que 
estão inseridas, e a criação de condições internas que 
lhes dêem ligeireza nas movimentações e exibilidade e 
qualidade na resposta aos desafios que constantemente 
lhe são colocados 

Para se desenvolverem com sucesso, as empresas têm 
de cuidar estratégicamente de três áreas: 

a) Tecnologia - As organizações industriais têm de 

analisar criticamente a tecnologia que usam em 

comparação com a disponível no mercado e a usada 
pelos seus concorrentes actuais é potenciais, de forma 

a tomarem decisões sobre os seus sistemas produtivos. 

Estrutura / Organização - Gerir a tecnologia e 

responder às solicitações do mercado em mutação, 

implica uma estrutura flexível, com diminuição de 
níveis hierárquicos e com descentralização de decisão, 
de uma forma coordenada e integrada. 

c) Recursos Humanos - tudo começa e acaba na 
qualidade dos recursos humanos disponíveis. Com a 
modificação das condições do mercado e da tecnologia, 
com esta a permitir a substituição cada vez em maior 
escala do homem pela máquina "inteligente" na 
execução de tarefas repetitivas e no seu controlo, são 
necessários recursos humanos com formação e 
capacidades distintas das exigidas num passado 
recente. Hoje, os recursos humanos têm de ter 
conhecimentos que lhes permitam, não só 
compreender e aceitar a mudança, mas também 
tomar decisões em tempo oportuno. 


Ni 


b 


A diminuição de níveis na estrutura organizacional, 
as consequências económico-financeiras das opções 
feitas na gestão dos sistemas produtivos e a utilização 
generalizada de meios computacionais nas organizações 
(tanto no sistema produtivo como no sistema de 
informação, gestão e controlo) têm, como resultado, uma 
cada vez maior interpenetração das componentes 
tecnológica e de gestão, nas tomadas de decisão. 


As organizações industriais precisam de téenicos que 
tenham não só formação tecnológica (conhecendo as 
tecnologias existentes e a sua evolução previsível) mas 
também formação em economia, gestão e ciências sociais 
(conhecendo o funcionamento do mercado, avaliando 0 
impacto de tecnologias alternativas, gerindo a 
tecnologia disponível, os meios humanos e materiais). 

A necessidade de um novo Lipo de quadros técnicos 
sentida em todos os países industrializados, no caso de 
Portugal assume especial relevância pela 
particularidade de algumas situações: 


- à integração na Comunidade Europeia e a progressiva 
abertura económica ao Mercado Único Europeu, 
trazem consigo, óbviamente, novas oportunidades 
mas, também novos perigos - um mercado muito mais 
alargado, uma competição acrescida - pelo que se 
torna indispensável que as empresas industriais 

tenham planos de desenvolvimento estratégico e 

formas de gestão dos meios disponíveis que lhes 

assegurem a sobrevivência e o desenvolvimento; 

na estrutura do tecido industrial são predominantes 

as PME's, com uma capacidade tecnológica e de gestão 

a necessitar de urgente reconversão e requalificação 

profunda (seja qual for o sector) ; 

- a enorme falta de quadros para gerir as organizações 
industriais de uma forma global e que integrem os 
conhecimentos da tecnologia e de gestão, que tenham 
a capacidade de se adaptar à mudança e os 
conhecimentos técnicos e atitude menta! de abertura 
que permitam ultrapassar o amadorismo e a rotina de 
execução. 


3 E as Universidades? 

As Universidades, e as suas Escolas, não podem 
portanto ficar indiferentes a esta revolução. Elas serão 
sempre, e estão a sê-lo, por elas afectadas. De uma forma 
mais evidente, aquelas que lidam mais de perto com o 
“mundo real”. 

Nestas estão claramente incluidas as Escolas de 
Engenharia. 

Assim, a questão que se põe também às Universidades 
é a de perceberem o que está a acontecer no mundo e 
adaptarem-se a essa mudança, de forma a ultrapassar os 
perigos e aproveitando as oportunidades. 

A Universidade, pela sua vocação de investigação e de 
formação de quadros, deve, não só adaptar-se à mudança 
como ter um papel importante de antecipação e 
condução. 

Mas para isso é fundamental! que, no que diz respeito 
à formação de quadros superiores, saiba quais as 
características do mercado onde estes irão funcionar e 
muitas vezes influenciar o seu desenvolvimento, 

Em particular, há que perceber o que se está a passar 
e se irá passar com as organizações industriais, e que 
tipo de quadros estas necessitam para vencer os desafios 
que a nova realidade lhes coloca. 

No que respeita a licenciaturas e tendo em atenção as 
mudanças no ambiente e a visão sistémica do mundo, 
devem ser fornecidos aos alunos conhecimentos que lhes 
permitam “compreender” esse mundo e uma formação 
técnica que lhes permita não só uma actuação eficaz no 
mundo do “trabalho” quando saírem da Escola mas 
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também a formação básica que lhes irá permitir uma 
adaptação constante. 
Assim devem ser contempladas na sua formação: 


a) Um sólido domínio das Ciências básicas 
"tradicionais" (matemática, física, química, etc). É 
indispensável que esta formação seja ministrada, 
pois, servirá de base e de “referência” à aquisição dos 
conhecimentos futuros, e permitirá uma maior 
capacidade de integração da formação continua do 
engenheiro. 

b) Informática. Os conhecimentos informáticos são 

fundamentais não só como suporte para O 

desenvolvimento eficiente das diferentes matérias 

mas, também, para que os alunos adquiram novas 
atitudes que serão de grande impacto na necessária 
transformação das organizações onde vão desenvolver 
as suas actividades profissionais, bem como na sua 

adaptação às transformações sociais a que estamos a 

assistir. 

Ciências aplicadas na area da sua 

“especialização”, e que corresponde à formação 

específica de cada “Lipo” de profissional. 

d) Formação na área de Economia, Gestão e 
Ciências Sociais. Não é possivel desconhecer duas 
realidades importantes: uma a de que são 
indispensáveis, no desempenho de funções técnicas, 
conhecimentos das áreas das Ciências Sociais e de 
Gestão com especial relevo para Gestão Global, 
Economia e Sociologia das Organizações, outra a de 
que, enquanto membros da sociedade, têm também 
um papel importante a desempenhar, pelo que 


Mi 


C 


aqueles conhecimentos são igualmente 
indispensáveis. 
4 Diferente realidade... o mesmo perfil 


de engenheiro ? 

Hoje, ao engenheiro não basta dominar a tecnologia é 
as técnicas operativas (sub-sistema tecnológico), é-lhe 
igualmente fundamental ter conhecimentos relativos 
aos sub-sistemas psicosocial, estrutural e de informação 
e gestão, sem o que não pode desempenhar, da forma 
mais eficaz, o seu papel nas organizações industriais. 

A actuação profissional do engenheiro ao ser 
enquadrado no actual ambiente das organizações não se 
limitará especificamente aos aspectos técnicos, mas, 
compreenderá também uma intervenção importante na 
área da gestão global. O engenheiro deverá ser um 
agente da indispensável mudança das empresas. 

Assim, o engenheiro precisa de uma formação que 
integre conhecimento das ciências básicas, de 
informática, de ciências aplicadas na sua área de 
“especialização” e nas áreas da economia, gestão e 
ciências sociais. 

Este novo perfil do engenheiro há muito que foi 
entendido internacionalmente como comprovam as 
seguintes citações: 


“Undergraduate education in engineering at MIT 
should prepare its graduates for leadership in 
technology, for professional excellence, and for rich 
lives of learning and reflection, through education in 


Luís Tadeu Almeida 


science and engineering with an emphasis on 

fundamentals, in essential partnership with the soctal 

sciences and the humanities, for the advancement of 
engineering and the betterment of society.” 

MIT, Goals for Undergraduate 

Engineering Education 


“Phe manager is the dynamic, life-giving element in 
every business. Without his leadership the "resources 
of production” remain resources and never become 
production. The quality and perfomance of its 
managers is the only effective advantage an enterprise 
ina competitive economy can have." 
Peter Drucker in 
"The Practice of Management” 


No IST, nas licenciaturas “tradicionais”, as 
componentes de Economia, Gestão e Ciências Sociais são 
diminutas ou inexistentes. Há assim, de modo a permitir 
satisfazer as necessidades apontadas, que introduzir 
alterações curriculares, que poderão passar por 
modificar os perfis/ramos existentes ou/e pelo 
aparecimento de novos perfis/ramos nas licenciaturas 
existentes, ou ainda pela flexibilização curricular 
através da instituição dum verdadeiro sistema de 
créditos, permitindo a "construção" de um currículo 
adaptável às necessidades e à vontade dos licenciados, 
sujeito a regras que assegurem uma sólida formação 
universitária, dentro do âmbito do IST, que 
naturalmente é também algo em constante adaptação. 

Na falta desta estrutura adoptou o IST uma política 
de lançamento de novas licenciaturas. 

Dentro desta política, começou a funcionar em 
1990/91 o 1º ano da Licenciatura em Engenharia e 
Gestão Industrial, que pretende ser uma resposta às 
novas necessidades de formação na área de gestão de 
tecnologia que se referiu anteriormente, de modo a 
contribuir para o ultrapassar das deficiências 
estruturais do tecido industrial português, formando 
engenheiros com uma visão interdisciplinar (sistémica) 
dos problemas industriais, e com sólidos conhecimentos 
de ciências básicas, das tecnologias fundamentais, de 
economia, gestão e ciências sociais, 

Esta licenciatura é assim um primeiro contributo 
importante do IST para vencer os desafios colocados pela 
Revolução Pós-Industrial, o qual não podia ser feito por 
outra entidade que não uma escola de engenharia, pois 
só nestas existem as condições básicas (o conhecimento 
da tecnologia) para que a integração das três 
componentes possa ser feita eficaz e mais rapidamente. 


õ As responsabilidades do IST face aos 
novos desafios 

O IST, sendo a maior Escola de Engenharia do Pais, 
tem por isso responsabilidades acrescidas, e o tipo de 
respostas que encontrar, para vencer os novos desafios, 
têm relevante importância no contexto de 
funcionamento e necessidades das organizações 
industriais. 

Um contributo importante foi dado pela criação da 
LEGI e pelo desenvolvimento da Área de Economia e 
Gestão. 
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Mas, para além da responsabilidade da LEGI, a área 
da gestão do IST deve ser um núcleo de ensino, 
investigação e prestação de serviços que, partindo do 
conhecimento das tecnologias disponíveis, permita uma 
melhor compreensão dos problemas enfrentados pelas 
Organizações Industriais, no que respeita à sua 
organização e gestão, bem como às formas de os 
ultrapassar. 

Neste contexto, à área de gestão cabe désempenhar 
uma dupla função: por um lado, incentivar nos alunos o 
desenvolvimento da atitude sistémica e por outro lado, 
fornecer-lhes um conjunto de instrumentos de análise e 
de resolução de problemas da Gestão das Organizações 
Industriais, que lhes permitam melhor compreender o 
mundo em que vão desenvolver a sua actividade 
profissional e torná-los mais apetrechados para resolver 


as siluações que ai vão encontrar, englobando as várias 
perspectivas em presença. 

Tem-se no entanto em conta que o objectivo do IST 
não é a preparação de técnicos em Gestão, mas sim 
engenheiros com uma visão de algumas questões de 
gestão com que vão ser confrontados. 

Mas este facto põe duas questões fundamentais 
interrelacionadas: 


- qual o papel da gestão numa Escola de Engenharia? 
- qual o futuro duma Escola de Engenharia como o IST? 


Quanto à primeira questão algumas respostas já 
foram dadas, e um processo de implementação está em 
marcha. 

Quanto à segunda, é necessário que as respostas 
sejam dadas rapidamente. A REVOLUÇÃO ESTÁ EM 
CURSO. 
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RESUMO 


Discute-se o conceito de "carácter iónico" e 
propõem-se métodos de o calcular, que se mostra serem 
mais precisos que o que recorre à tabela de 
electronegatividades de Pauling. 


The concept of "ionic character” has been differently 
defined, and as a result, it has been misleading to reason 
in terms of it, unless one really explicits the precise 
context in which the "ionic character" is discussed. 


(A) Pauling [1] has suggested that the experimental 
dipole moment, u, of a heteronuclear diatomic 
molecule may be used to estimate the ionic character 
of a single bond. Viewed as two point charges (+q) 
and (-q), separated by the internuclear distance, d 

u=Iglxd 
where Igl =6 Xe -— e=4.803 XxX 1010esu. 
So, for any alkali metal halide molecule, the ionic 
character is given by (ô X 100), which corresponds to 
the fraction (expressed in percentage) of the electron 
charge present in each electric pole. 


(B 


Di 


The above definition of ionic character implies that a 
purely covalent bond would have no dipole moment. 
However, owing to a difference in shape or size of the 
two overlapping orbitals, some distortion might give 
rise to a dipole moment even though the bond pair is 
equally distributed. Coulson [2] conceives the ionic 
character as 10042/(1 + A2) where 
P = Peov + À Pion; 

W is the complete molecular wave function, and À is 
the Coulson constant sometimes called the degree of 
ionicity. 


(C) Since Pauling [1] postulated that the 
electronegativity difference, for a single covalent 
bond between two atoms X and Y, measures the 
extent to which the bond in this XY molecule is 
polar or partially ionic in character, it has been 
accepted that the partition of the single bond 
enthalpy, 4H(X-Y), leads to a covalent term and a 
ionic term. 

First, the covalent bond enthalpy term, 4Hcou X-Y), 
was identified with the arithmetic mean of the 
single bond enthalpies of the molecules Xy and Ya. 
Later, it was suggested that 4H .gu(X-Y) would be 


Original submetido para publicação em 19 de Agosto de 1957 


ABSTRACT 


The concept of “"ionic character" is discussed. 
Methods for its evaluation are proposed. They are 
shown to be more precise than the estimation of ionic 
character through the Table of Electronegativities of 
Pauline. 


better estimed through the geometric mean. 
Obviously one could infer that 

AH ionlA-T)=AH(Ã-Y)-AH ou AY). 
Let us work a particular set of exampies (potassium 
halide molecules) that happen to be electrolyte 
diatomic molecules. + 


Diatomic Molecule (From refs. [3] and [4]) 


E TESE 
| (kJ mol-?) 
MRE RC 
KAMEN ERR 
EEE ESEC 
923 | 


lonic Character 


Estimated from 
Table of Electro- 
negativities 
(Pauling/Allred) 


po 


Diatomic 


u x 100 
Molecule 


AH ud X-Yix 100 


4.803 xd AHX-Y) 
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(D) A different model of deriving the ionic character of a 
heteronuclear molecule is inspired in the extended- 
calculations of the lattice energy of alkali metal 
halide crystals, that are based on the Coulomb, 
Born, London and zero-point energy contributions 


õ 
- 
x 
o 
4 


[arde fram [mo [rs [res 
ie | to (to | tio at 


713 


p = internal energy percentage of ionic nature 


Since the experimental value of the lattice energy 
for any alkali metal halide crystal is higher than Lhe 
extended-calculated value, this suggests that there 
is a systematic error in such estimations due to 
disregarding a covalent term. 


For diatomic molecules we may apply similar 
equations to Born-Landé and London equations 


[5 xr, (1) 
AU g+c KV 9) =—1xa XAU se XY) 
I à Bam 
ra AD TE AE ) 
(4ne,) dº 2 (E +E O) 
where 


A=Madelung constant for NaC] 
structure = 17476 

a + = polarizability of the cation 

= (0.93x 10-40 F m2 
-= polarizability of the anion 

E + = ionization energy 0[K+ = 3.051 kJ mol-l 

E = modulus of the electron affinity of the 
correspondent halide atom 


type of 


av, 


a 


It is obvious that the ionic contribution to cohesive 
energies is greater in the solid state than in the 
molecular gaseous state. 


Moreover the energy percentage of ionic nature of 
the molecules is closest to the enthalpy approach of 
ionic character of the same molecules, 


So, Summarizing, 


- one could notice that the concept of "ionic character” 
may be understood in terms of charge distribution in 
the molecular orbital y (Coulson), dipole moment 
(Pauling) or bond energy; 


- the bond energy ionic character may be estimated 
either through the Pauling equation or by extended- 
calculated contributions; 


- the estimation of ionic character through the Table of 
electronegativities of Pauling/Allred is very rough. 
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RESUMO 

Apresenta-se e discute-se teoricamente uma técnica 
não convencional de controlo baseada em regras 
"fuzzy". E feita à sua aplicação a um problema 
concreto de engenharia. 


| Complexity, Fuzziness and Fuzzy 
Control 

Since the 50's that a great deal of the contemporary 
societies, especially the Western ones, have been 
undergóing the influence of the a Renaissantist-like 
movement (that many claim to be the ushering in a new 
information age”) which, arising almost simultaneously 
in various areas of knowledge (particulary in those 
where control, information, computer technology and 
related fields have been especially concerned) has, in 
turn, begotten the emergence and development of new 
fields of inquiry such as systems science, cybernetics, 
information science, artificial intelligence, etc. 

Al of these changes led to a profound reappraisal not 
only of the entire paradigm on which the whole scientific 
panorama had, until the 50's, been based but also to the 
discover and/or reformulation and extension of many of 
its researched tenets. From the manifold, albeit cognate, 
subject matters which the emergence of these new 
themes entailed there are three which, due to their 
present-day importance, deserve a rather especial 
attention. 

The first of these subject deals with complexity., 
Advances in computer technology have been extending 
our capabilities of coping with systems of an 
increasingly broad range, this extension including 
systems which due to their nature and complexity have, 
until quite recently, been regarded as intractable. 
However, insofar as the level of complexity we can 
manage continues to increase, also we begin to realize 
that there are some fundamental limits which, by no 
means at all, can be surmounted. Or, putting the matter 
differently, we also begin to understand that the need for 
simplification is often unavoidable. One of the processes 
used for simplifying complex systems - the most 
significant one, perhaps - is to allow some degree of 


This paper is supported by the Nato Special Programme Pannel on 
Advanced Educational Technology and by CAUL (Centro de 
Automática da Universidade de Lisboa). None the less, Lhe 
responsability for the views and opinions contained herein pertains, 
entirely, to the authors. 


Original submetido para publicação em 18 de Junho de 1991 


ABSTRACT 


A non-conventional control technique based on 
fuzzy rules is presented and theoretically discussed. 
The analisys of the case study illustrates its 
application to a concrete control problem. 


uncertainty in its description. Generally speaking, a 
good simplification must minimize the loss of 
information (regarded as) relevant to the problem which 
is being concerned. And because this loss implies 
(amongst other things) an increasing degree of 
uncertainty, then it follows that complexity, 
information and uncertainty become matters mutually 
related to. 

Later we shall return to this assertion. For the 
moment, taking 'relevancy' as a starting point (How can 
one know what is 'relevant' or not for the solution of 
some problem?), our attention will be directed to the 
second of the subject matters previously referred to. It 
deals now with knowledge. 

Together with such advances in computer technology 
this is, perhaps, the area which more transformations 
has undergone in the last fourty years. As a matter of 
fact, from a period in which emphasis was laid on an 
attempt to know the character of the 'external" (or 
Physical") world we have now reached a different (or 
reflective) stage in which the main focus of interest is, 
rather, devoted to what may be called 'the knowledge of 
the knowledge". In other words, we want to know not 
only 'facts' or 'truths' as before but also what we can or 
cannot do with this knowledge, what can or cannot be 
known, how we know and how to cope with knowledge 
(and ignorance) itself. One of the major consequences of 
this shift (closely related to the methodology which the 
use of computers entails) is a renewed amount ofinterest 
on psychological phenomena. This interest has been 
translated in the emergence of new approaches to mind's 
functioning, to its architecture, to reasoning, to new 
logics, etc.; briefly, to new trends in the psychological 
sciences, of which the so-called 'cognitive psychology" is 
the present-day utmost representative. 

We shall not dispute the fruitfulness of this cognitive 
current in the understanding (and further simulation, 
replication even) of some mental phenomena, 
particularly of those which Jackendorff asserts to 
pertain to what he calls "the computational mind", But, 
in our opinion,such an overall approach tells us only a 
part of the story. Human beings are far more complex 
and puzzling than simple serial or parallel logic 
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machines, the word “logic” being here referred to its 
traditional meaning. That we have minds and brains no 
one denies. But we do have bodies too and the model of 
the brain as a mass of neurons-like-on/off electronic 
switches controlling the body not only has proved to be 
erroneous in the light of the most recent advances in 
neurobiology but has also left unsolved the relationships 
between brains and minds as well as between minds and 
bodies. As a further point, the emphasis laid on rational 
and cognitive processes alone so characteristic of that 
approach (excluding affectivity and irrationality almost 
by edict) introduces a false and rather artificial division 
into the mind's functioning, far from what is presently 
know about normal and abnormal mental processes. 
Even the modelling of (classical) logic reasoning - one of 
the strongest bastions of mind-like simulations - 
undergoes severe criticisms if compared with Lhe real 
processes we use to get some 'neat' conclusion from 
previously known (and also supposedly 'neat') premises. 
As a matter of fact, beyond the ubiquitous and, yet, 
obscure problems which the (selfjorganizing process that 
the coherence of our thoughts/utterances beset, a great 
deal of the terms/concepts/inference rules we use 
whenever we are engaged in conversation either with 
other people or even with ourselves is far from being 
'neat' and/or 'objective' (as the traditional logicians 
claim). On the contrary, and to a great extent, they are 
essentially vague, they usually comprise several types of 
(Lime-variable) self-referentials and, consequently, they 
often entail a large amount of subjective and imprecise 
appraisals in which the individuals 'experiences (as well 
as their idiosincrasies) play a non-neglectful role. 

This vagueness, subjectivity and one's individuality 
underly, precisely, the last of the three subject matters 
previously referred to. For the sake of conciseness we 
shall name it the gap between hard' and 'soft' sciences, 
this expression entailing, broadly speaking, Lhe gap 
between the high standards of precision which, on the 
one hand, prevail in mathematics, physics, chemistry, 
engineering and so forth, and Lhe imprecision of which 
much of sociology, psychology, education, political 
science, history, philosophy, linguistics, anthropology, 
literature, art and related fields are, on the other hand, 
pervaded. There are manifold reasons justifying the 
sharp contrast between Lhese two Lypes of precision 
standards!. But, from them all, the deeply entrenched 
tradition of precise and quantitative reasoning in the 
“hard” disciplines (Lhe roots of which can be directly 
related to the effective success of the mechanistie 
paradigm in the descrition/explanation/prediction of a 
great many deal of relatively non-complex "natural" 
systems) plays, undoubtedly, the most prominent role. 
As a matter of fact, this (possibility of) quantification 
establishes a clear difference as regards Lhe "soft" 
sciences, whose main bias is centered on systems where 
those qualitative factors intrinsically associated with the 
complexity of the human functioning (this including 
human values, norms, judgements, perceptions, 
emotions and so forth as well as human interactions) 
cannot, by no means at all, be disavowed or neglected. 


1 Seel13 |for more details.. 


Due to the enormous influence which the mechanistic 
or “classical” paradigm is still exerting upon the 
present-day scientific panorama, the conventional 
mathematical techniques have been, and continue yet, 
to be applied to the analysis of such humanistLic systems. 
In spite of this, it has also becoming increasingly clear 
that the complexity, particular characteristics and 
specific problems which such systems offer, call for 
approaches that have to be significantly different both in 
spirit and in substance from the traditional methods. 
One of the approaches which, in the last two decades, 
has progressively been recognized as providing solutions 
for a great deal of such problems is the so-called theory of 
fuzzy subsets of Lofti Zadeh. 

The theory arose explicitly in the 60's 2, However the 
french expression "ensemble flou" had already been used 
in 1951 by Menger when, based upon a probabilistic 
interpretation, he set up whal is presently named the 
max-product operation of a fuzzy transitive relation. 
From a semantic perspective, Zadeh's approach is also 
rather close from Black's work on “vagueness"” (Black, 
1937); his "consistency profiles" can, indeed, be regarded 
as the ancestors of the present-day "fuzzy membership 
functions” which, since the inception of Zadeh's theory, 
are used to describe words, concepts or situations which 
the classical approach is unable to deal with. 

Some of these situations recover, precisely, the three 
major subjects previously referred to, 

As regards complexity and knowledge, because it has 
been recently realisedº that there are several 
fundamental types of uncertainty (which, together with 
vagueness, also beset various types of ambiguity such as 
non-specificity, dissonance and confusion) each of which 
plays a distinct role in the respective simplification 
problem. As a further point, the mathematical 
framework into which these uncertainties can be more 
properly inserted has also proved to be available in 
terms of an overall theory of fuzzy sets and fuzzy 
measures. 

As regards the debate hard/soft sciences because a 
great deal of the (potential) conflicts between them - 
particularly those referred either to the semantics of 
propositions (excluded from the neutral, syntatic 
language either of General Systems Theory or, more 
generally, from all the disciplines which obey the 
classical paradigm of physical/quasi-physical 


disciplines) or to reasoning based upon imprecise 


premises - have both been overcome by a fuzzy extension 
of Zadeh's theory to a new field of inquiry, the so-called 
“approximate reasoning". We shall not go into details 
about its tenets. For the moment we shall only stress out 
that this approach is based on two main concepts, a) that 
of linguistic variable, say of a variable whose values are 
words or sentences in a natural or artificial language 
(for instance, "age" may be regarded as a linguistic 
variable if its values are assumed to be "young", "not 
young”, “very young”, "old", "not very old”, etc) and b) 
that concerning the inference rules (nothing but 
particular fuzzy operations) which allow to get some 


2 See [17] 
3 See [7,8] 
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conclusion from premises which are neither very exact 
nor very inexact. 

Fuzzy sets and fuzzy logics underly too not only these 
themes but also the areas where computer simulations 
of the interactions either between human beings, 
between human and machines or even between 
machines, are concretely involved. One of these areas 
(where complexity, uncertainty, linguistic variables and 
linguistic rules as well as interactions Man/machines 
are altogether related) deals precisely with the major 
objective of this paper: the theory of the so-called 
linguistic controllers, an engineering control technique 
based upon fuzzy concepts and fuzzy operations which is 
becoming increasingly used in present times. 

The first concrete application of Zadeh's fuzzy theory 
to control problems is due to Mamdami and Assilian 
(Mamdami & Assilian, 1975). In a classical (and now 
historical) paper entitled "An Experiment in Linguistic 
Synthesis with Fuzzy Logic Controller" they 
demonstrated how the main theoretical fuzzy operations 
can be used to the effective modelling and control of 
systems which in that epoch were incapable of being 
adequately treated according to the conventional 
techniques. Since then, Mamdami's revolutionary 
approach has given rise to more Lhan 3000 papers on 
fuzzy control applied to several, distinct areas. 

Despite this enormous amount of work a great deal of 
the present-day literature about the theme undergoes, 
from the authors! perspective, some severe criticisms one 
of which is especially poignant. As a matter of fact, a 
survey of such works (until 1990 at least) shows that 
there is an enormous gap between the type of papers 
where extremely specialized control problems are the 
main focus of attention (the bias of which as well as their 
understanding is only restricted to some few experts in 
such specialized areas) and some other papers/books 
which, albeit related to fuzzy control, deal with so 
general notions that their application to a control 
question (out of the specialized publications) becomes, 
for the non-experts, extremely difficult. The overcoming 
of this gap is one of the main challenges of the present 
work. 

In order to solve this (and other) problem(s) we shall 
use in the following considerations an approaching 
strategy which will be divided into two main phases. In 
the first, we shall briefly survey some of the major 
questions which the modelling and simulation of control 
systems according to the conventional perspective 
usually entail. The main structural and methodological 
differences existing between conventional and linguistic 
controllers will, next, be pointed out. In the second 
phasis, a great deal of attention will be devoted to the 
various stages which the design of any linguistic 
controller (whatever it may be) must follow. A concrete, 
albeit simply illustrative, example will provide the 
reader with a computational support for further 
developments. 


2 Conventional and Linguistic 
Controllers: Main Differences 
Since the very beginnings of control theory that the 
designers of automatic control systems have had to face 


some typical types of problems, two of which may be 
characterized as follows. In the first, the designers' 
objective seeks to hold the value(s) of some “output” 
variable(s) within a certain predefined and 
predetermined region, this taking into account a given 
set of possible disturbances against which Lhe 
(endeavoured) control system must react. In the second, 
that objective aims to pursuit (or "follow") a certain, 
usually moving, target with the maximal accuracy so 
that the actions to be performed by the control system 
can both antecipate and correct some (possible) set of 
“unpredictable” target's movements. 

According to the conventional perspective both of 
these control systems obey Lhe general scheme depicted 
in fig.1.We shall not go into details either about the 
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fig. 1 

particular forms that this scheme may take in terms of 
concrete applications or about the specific mathematical 
techniques they entail, attention has, however, Lo be 
paid to its underlying philosophy. As a matter of fact, 
the controller's design for a specific application requires 
three (at least) stages which, respectively, deal with a) 
the analysis of the actual system, b) a decision criterion 
based on the linearized model, on some desired state and 
on some error criterion function, c) the synthesis of the 
controller's struture as well as of its parameters. lt 
follows, consequently, that such conventional 
controllers are able to provide accurate control 
whenever the actual systems are (relatively) 
non-complex, they can be matched to linearized models, 
the desired state can be well-defined, the error criterion 
function is easy enough to cover all the system's 
purposes and the quality/amount of the measure of the 
parameters' values is high. Í 

A totally different situation occurs, however, when 
the process or system to be modelled is extremely 
complex (leading, for example, to deep uncertainties 
about its internal structure), presents non-linear 
relations between its variables/parameters, these 
parameters are time-variable (sometimes with 
ill-defined time-variations), the system contains 
non-specified time-delays and multi-objects to be 
controlled are effectively needed. 

Two possible ways of overcoming these problems have 
been used. Firstly, the designer constructs simpler and 
simpler replica of the system (according to some 
predefined set of simplificative hypotheses) using for it 
(more or less) complex structured algorithms. Since, 
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however, these algorithms translate Lhe knowledge 
which the designer has of the system and this, by inicial 
assumption, is extremely complex (which implies a 
certain degree of uncertainty/ignorance from the 
designer's side) then, whenever a non-predicted 
alteration of some var'able/parameter/subsystem of the 
global system occurs or its structure changes, it follows 
that new and new reconstructions of these algorithms 
are always mandatorily required. Which, in turn, will 
imply an additional (and often tremendous) amount of 
work whenever these situations arise. 

The second way of cvercoming such problems takes 
into account the natur:! lexibility, adaptativeness and 
expertness of humar cperators. In this case, however, 
the controller's designer has not only to provide the 
operator with a set of “instructions” but to consider also 
that this operator is able to interpret and understand 
them in a language which, because founded on the 
natural language, contains terms which are highly 
vague and imprecise. 

To illustrate this point let us consider a concrete 


example. fig 2 shows an extract from American textbook 
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| action b. 


«Increase fuel To maintain burning 


rate zone temperature 


« Decrease fuel | To raise percentage DÊ 


rate slightly DEygen 


- Reduce fuel To increase percentage 
rate 
-« Reduce I.,D. 


Fan speed 


To lower bach-end teme 
perature and maintain 
burning zone temperatu- 
re 


Increase I.D.| To raise back-end tem=- 


Fan speed perature 
Increase fuel |To maintain burning 


rate 


(zone temperature 


HOHE. 


not get avergon- 


However do 
Fident and keep 


ali conditions 
under observation 


tig.2 


for cement kiln operators*. 
This extract describes the actions Lhat must be 
executed by a human operator in order to adjust the fuel 


4 See [5] 


of oxygen for action b. | 


rate, kiln speed and air volume through the kiln under 
several variable conditions characterized by the 
temperature in the burning zone (BZ), oxigen 
percentages in the exhaust gas (0X) and the 
temperature at the back-end of the kiln (BE). 

Beyond the actions which the operator has to perform 
(meaningless in this paper), the aspect here which must 
be pointed out lays emphasis on such terms as "low", 
“high”, "OK", "slightly increase”, "slightly decrease” 
and so forth which, despite their use in our everyday 
utterances, have a content which, by no means, can be 
translated according to the “exact” representations of 
the conventional techniques. They can, none the less, be 
expressed in terms of fuzzy proposttions and fuzzy rules. 
And because these expressions have, since Mamdami's 
paper, been able of computer translation then, it follows 
that fuzzy or linguistic controllers are a good means 
either to substitute human operators in control 
situations so far regarded out of court or to solve control 
problems where complexity, non-linearities, etc. are 
prevailling. 

This being assumed the general scheme of a linguistic 
controller (depicted in fig.3) can now be easily 
undeerstood. 


FUZZY INFERENCE ENGINE 


Pressanicescnanemas 


| CONTROL 


5 og 


chi Rildo X, ES. + AND Hg 1a o + THEN WITH. . 
| Ê : 
| PURPDBES + 
| CONTROL ACTUAL 
DES IBRED 
“ERROR” COMMAND STATE 
STATE 
bx Y x 


polias Ep 


DECISION CONTROLLED | 
UNIT SYSTEM 


fig. 3 


The system contains a fuzzy control unit where a set 
of control rules R, is implemented. These rules may be 
described in the simplest case as "If x is A; then y is By" 
where x corresponds to a linguistic evaluation of the 
system's state and y (the control command) is the action 
the system has decided to execute in order to satisfy the 
controller's objective. Since this objective coincides with 
the designer's purposes then it can also be asserted that 
the control unit materializes the designer's reasoning in 
a way which, because vague and imprecise (albeit more 
adaptative and flexible) is much closer of the human 
functioning. 

These (and other) features of fuzzy controllers 
establish a clear cut as regards the conventional contro! 
techniques, particularly those in which "expert systems" 
(in the sense of Al) are used. As a matter of fact, despite 
their (apparent) similarities in terms of decision and 
implementation of rules, the remaining methodologica! 
and even philosophical differences are, however, so 
poignant that only a detailed description of how u 
linguistic controller works is powerful enough to 
account for them. Such a description (nothing but the 
proper design of the controller) will be the content of the 
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next considerations. 


3 Design of a Linguistic Controller: 
Main Stages 

Let us assume that the designer has some previous 
knowledge of how the process/system to be modelled 
broadly works. In other words, let us suppose that he/she 
has already made a distinction between some of the 
input and output variables and that he/she knows their 
respective numerical values, particularly the minimal 
and maximal values that these variables can take 
throughout an enlarged observation period. Based on 
this previous knowledge the design and effective 
implementation of a fuzzy controller relating such 
(selected) input and output variables obeys six main 
stages which, respectively, deal with: 


-the definition of the universes of discourse 
corresponding to the (hypothesized) input and output 
variables; 

- the fuzzification of these universes using fuzzy terms 
("great", "small" and so forth); 

-the setting up of the partial 
relationships (fuzzy inference rules), 

-the construction of the global input/output 
relationship; 

- the defuzzification stage; 

- the analysis of the obtained results (which may lead 
to the introduction of new universes of discourse, new 
inputs, new rules, etc.) until that such results are 
assumed to satisfy the designer's objectives. 


input/output 


A more detailed description of these stages (as well as 
of some of the problems they entail) is provided in the 
sequel, 


3.1  Universes of discourse 

The knowledge of the universe of discourse of every 
(hypothesized)variable is a crucial piece of information 
for the construction of the controller. Manifold reasons 
justify such an importance but they will be more clearly 
understood as far as the remainder phases of the design 
will be explicated. For the moment it suffices to bear in 
mind that a change in the range of these universes will 
influence pratically all the further implementation 
stages (this including the definition of curves which 
translate the terms used to characterize the controller, 
the inference rules, the mathematical operations they 
entail, the type and number of the discretization stage 
that such operations require, etc.). 

As previously emphasized the boundaries of these 
universes are defined taking into account the minimal 
and maximal values which the (selected) input and 
output variables may assume throughout a more or less 
extended period of observation. These boundaries are 
usually, normalized (tipically between O and +11 or, 
sometimes, between -1 and + 1). Appropriate scaling 
factors (not necessarily alike for all the variables, their 
universes of discourse or, if this is the case, for the 
positive and negative parts of such universes) must, 
therefore, be chosen. Because irrelevant for our present 
purposes these scaling factors will, however, be excluded 


from these theoretical considerations; thus, only 
normalized universes will, thenceforth, be considered. 


3.2  Fuzzyfication of the inpuí(s) and 
output(s) 

Once the universes of discourse of the (hypothesized) 
variables are defined the designer must next translate 
the observed variation of their numerical values in 
linguistic terms. Which, putting the matter differently, 


LOW 


MEDIUM HIGH 


=) sind [y =*t Alf, at 
“44 a (x Ergsinção (x p) ) 
a<x<b, an<zx<hb, 
= 0, b <a. =), d<a. mel; bh = x. 
where: na= 0.1, b=05 a=0.2, b=0.45, c=0.55, d=0,6 a= 0.5, b= 0.9 
a) 


12 
& (0.5) 150 


uix)j= e k>0 uix)=e pix) = 0, U<ar<a, 
k=5 k= 5, ko, 
=y-ettrals es, 

b) 0=0.6, k=50 


ge — k >A, “piel. >], x) =0, pexca, 
iam Lx? d PTS ei me k(x-a)É = 
= Lk(xa)? m<x<i, 
E=3 k = 50 | 
€) o=06, k=2H0 
fig. 4 


is tantamount to saying that the designer must express 
his/her subjective "estimation" of such variations in 
natural words as “great”, "small", "slightly increase" 
and so forth. Since these terms (also named 
"qualificatives", "attributes" or "linguistic variables") 
are equivalently translated by appropiate fuzzy 
membership functions then the type and number of such 
“linguistic estimations" will thenceforth define not only 
the maximum number of rules that the control unit 
must contain (in principle, for a SISO system, this 
number equals the product of the number of input 
attributes by the number of output attributes) but also 
(once implemented) a first appraisal of the behaviour 
that the system must futurely exhibit, 

The curves by means of which such membership 
functions are graphically displayed can take several 
forms. fig.4 depicts some of the most typical ones. 

It is worth noting that all these curves depend on 
parameters ("k”" or "a" and "b") whose values may 
undergo several variations either according to the 
designer's perspectives (for instance, as far as his/her 
knowledge of the system's functioning increases) or, in 
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more elaborate controllers where "adaptation" and 
“learning” have (automatically) been introduced, taking 
into account some optimization criterion. 


3.3 Construction of the fuzzy engine 

Ônce the linguistic values of the input and output 
variables (as well as of their corresponding curves) have 
been selected, the next step of the controller's implemen- 
tation consists in the setting up of the inference rules 
which superordinate the relationships between input 
qualificatives and desired output attributes. 

Two different, albeit cognate, types of problems frame 
this stage of the controller's design: one deals with all 
the mathematical paraphernalia that such relationships 
beset, the other the questions which the (necessary) 
discretization of the already selected universes of 
discourse may eventually beget. Both of these problems 
will be more circunstantially approached in the sequel. 


3.3.1 Partial fuzzy rules 
Let P;(i=1,2,3,..n)and Q4=1,2,3, ...m) be fuzzy 
propositions including fuzzy predicates such as 


Pr="XisAto Qu="risBro MM 
L ! J J 


where X and Y are linguistic variables and A, and B, 
represent their values or attributional properties (such 
as "great", "medium", "small", "low", "very low", "high", 
“moderately high”, etc) assigned to X and Y according to 
the designer's criterion, Then, for a certain value of i and 
j, say k, a linguistic rule is a meta-implication 
statementô which, in the simplest form, is expressed by 
“IF P, THEN Qk' by means of some implicit relation Ryp 
(which is said to matertalize the rule) 


Rk 
Pr; — Qk 
which, in a more suggestive graphical way may depicted 
as 


P, Qy 


ami 


fig. 5 


“ The reader must carefully bear in mind the distinction between 
implication and meta-implication. Às à matter of fact, let P and Q 
be two formal propositions and let us denote implication and 
meta-implication respectively by the symbols ">" and ">". Then 
the composite proposition "P implies Q" corresponds to the truth 
table depicted bellow;"F" and “Tº being the usual “values” of 
“falsehood” and “truthfulness". As a further point, if À and B are 
sets respectively related to P and Q then such a composite 
proposition 15 given by the set theoretical operation AUB, On the 
contrary, in the composite proposition "P meta-implies Q" it is 
assumed that whenever P is true then Q is also supposed Lo be 
true, although nothing can be inferred from the meta-implication 
if P is false. The meta-implication introduces yet two additional 
features, a cause/effect relationship and the content of Lhe proper 
propositions, both of which are not entailed by the implication. 
For more details see [6]. 


P Q P>Q 
b b T 
F I T 
T F 3 
Tr 1 F 


In linguistic terms P, expresses some propositional 
condition to be satisfied (antecedent) in order that Qk 
(the consequent) can also be verified. Obviously, when 
translated in the controller's "language", Lhe 
'antecedent' corresponds to the conditions which the 
input is supposed to satisfy and the 'consequent' the 
conditions which the output has to verify from the 
designer's perspective. 

This being assumed the rule which materializes the 
relation Ry between Pp and Q, is usually obtained by 
means of Zadeh's inference operation 


o Tn dó 


cs T (2) 
=MINI(a,) 56, 


where a, and by are the (computational) vectors by 
means of which the (discretized) membership functions 
that translate the linguistic values A, and B, are 
expressed, (as)! means the transpose of a, and the 
operation "w"” symbolizes (in this case) the so-called 
direct (or cartesian) product. 

More complex rules entailing several, distinct, 
antecedents P4;, Pa; ...., and one consequent Q, (1,J, ....= 
| Za mM sTE Il, 2 «8) logd to the composite 
meta-implication statement " IF Pj, AND Ps; AND 
EE - THEN Q,". Again, for some specific k, the implicit 
relation Ry may be expressed by 


R, =P, ANDP, AND..)xQ,]= 


Es T (3) 
= MINI(a jta,,t..) ,b,] 


where ajk, a9k, ... corresponds, as before, to the specific 
(discretized) membership functions which linguistically 
describe the attributional properties Ay; €P1;, Ag;€ Po, 
etc., “x” is the inference operation previously defined 
and "t" is a triangular norm (or "t-norm”) which obeys 
the following properties: 


l)tisa relation [0,1]x[0,1] [0,1] (4) 
For any membership functionsa,b,c,d €[0,1]: 


2) t is a non-decreasing operation in its 
arguments, say, 


f(a<=b)and (d=c) 
then (atd)j< (btce) (5) 
3) tis a commutative operation, 
atb=bta (6) 
4) tis associative, 
(atb)tc = atíbtc) (7) 
9) tsatisfies the boundary conditions 
at0=0 (8) 
atl=1 (9) 


6) In all cases 


atb < MIN (ab) (10) 


Examples of these t-norms are provided by: 


ajatb = MIN (a,b); 
b)jatb = a.b where"."is the usual scalar 
multiplication; 
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c)atb=1-MIN[L((1-a)w + (b)w)lw] with wz=; 
d)Jatb= MAXI[O,(A+1)(a+b-1)-Aab] with 
A=-l; 

etc. 

From these t-norms the most used is, undoubtedly, a). 
However, c) which corresponds to the so-called Yager's 
intersection, may provide stronger or weaker 
intersections as far as w varies (a property which, 
broadly speaking, is shared by all the t-norms which 
depend on some parameter). This feature can be 
advantageously profited in the design of self-learning 
fuzzy controllers. 


3.3.2. Discretization of the universes of 
discourse 

The mathematical operations which the construction 
of each one of the preceding rules beset require, 
mandatorily, the discretization of the different 
membership functions by the means of which attributes 
are translated. À great deal of problems is encompassed 
in this discrezation. Two of them are here especially 
poignant, a) the number of intervals into which the 
universes of discourse of every relevant variable must be 
divided and b) their respective ranges. The first rises a 
question which takes into account both the accuracy of 
the endeavoured results (which requires a greater 
number of intervals) and the computational effort 
(which, in turn, demands a redution of such a number). 
A possible way of overcoming this ambiguity consists in 
the use of discrete intervals with different ranges: near 
the contolled region this range is made small, away from 
it the range is increased. 

Whatever may be Lhe selected option (variable from 
case to case) the designer has, however, to take always 
into account the computer's speed as regards the set of 
possible disturbance to which the whole system must 
react. Hence, the number of selected intervals, their 
ranges, the type and characteristics of such possible 
disturbances and the proper features of the computer 
used becomes matters mutually related to. Nonetheless, 
insofar as we know, no specific research study has been 
specifically performed upon this particular, albeit 
crucial, aspect of the controller's design. 

As a pratical advice it is usual to begin with a certain 
number of equalized intervals (satisfying some rough 
accuracy criterion) and next to improve the results 
obtained either by increasing that number or (if the 
controller region is already well-defined) by changing 
their ranges; this, taking always into account eventual 
delays between the controlling actions executed by Lhe 
fuzzy engine and the real disturbances it has to correct. 


3.4  Globalinput/output relationship 

Assuming that this problem concerning the 
discretization of the universes of discourse is solved and 
that all the partial rules have already been defined the 
next stage of the designer's strategy will require the 
setting up of the global relation Ry which entails the 
different R; (i=1,2,...N) established in Section 3.3.1. 

Broadly speaking Ry corresponds to an operation 
which join the R; according to two main sentences 


ALSO 


(F...THEN..), UF...THEN ..), 


asda dd 


ELSE, ELSE 
The interpretation assigned to these “also” or “else” 
leads to various types of operations. The most usual 
corresponds to the classical fuzzy union 


R.UR.=MAAIr.r.) 
i J a 


ER: (11) 
L Ê 
r ER 
J J 
None the less, and similarity with the extention of the 
classical fuzzy intersection by means of t-norms, it is 
also possible to define a generalized union or s-norm 
(also named t-conorm). These s-norms obey all the 
properties previously referred to t-norms except that one 
which deals with the boundary conditions which, now, 
verify the expressions 
(12) 


II 


rs 0=r 
Ê | 


(13) 


It can also be demonstrated that for any 
corresponding r;, r; 


MALA lr.r)Sr sr. 
Md b J 


Prgl=] 
| 


(14) 


As in the case of t-norms, various types of s-norms 
may be found in the present-day literature. Some typical] 
examples of these norms are given bellow for two 
general membership functions/membership relations 
"ar and “hr. 


a)asb = MAX (a,b) 
b)asb=1-MAXIO,(1-a)w+(1-b)w-1] with 
w € (-c0, 00) 

cjasb=a+b-ab 

d) sb= MIN[1L,a+b+Aab] with Ã=-1 

e)asb = MIN[1, (aw+bw)lw] with w É (0, 00) 

etc. 

Again we remark that values of w or À may lead to 
weaker or stronger values of Lhe generalized union, a 
property which may be advantageously profited in the 
construction of adaptative, self-learning controllers. 


3.5 The defuzzification stage 

Based on the knowledge both of Ry and of the 
linguistic attributes corresponding to some numerical 
input value, another fuzzy operation (max-min 
composition, etc.) will allow, next, the specification of 
the global linguistic value for the output variable 
corresponding to that input value. 

Tipically, as exemplified in fig.6, this global output is 
given by the composition (fuzzy union or some of the 
operations which s-norms encompass) of several 
“weighted" outputs, each of which translates the 
influence of the particular rules that such an input 
value entails. 

In order to obtain a definite numerical value for the 
controller's output, a defuzzification procedure must 
therefore be selected and implemented by the designer. 

There are several possible criteria ruling this 
selection. Three of them, the most commonly used, are 
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fig. 6 
described in the sequel, 


3.5.1 Max-procedure 

This procedure relies on choosing the output value as 
being that one for which the membership function 
attains a maximum. Hence, il yy is assumed to be this 
value and Y (y) represents the global output 
membership function, then we have 


MAX YO)=Y(yy) 


Although simple, the method undergoes two severe 
restrictions. On the one hand, it may beset ambiguities 
as in the case where more than one element of Y reaches 
that maximal value. A possible way of overcoming this 
situation consists in finding their average value 
(mean-maxima approach) which will futurely work as 
the desired output value. Hence, for n maximal values y'! 


(15) 


R y' 
= > (>) 
0º 4a (16) 
= 
On the other hand, the procedure does not also 
account for the shape of the output membership 
functions, an information which, in some situations, 
may be crucial for the understanding/correction of the 
further control actions. 


3.5.2 Centre of gravity procedure 
In this case the value yp is determined according to 
the usual estimating procedure in which the 'center of 


gravity' of the whole surface which Y(y) encloses must 
be found. Thus 


| yY (y) dy 
Y (17) 


| Y (y) dy 
r 


3.5.3 Indexed centre of gravity procedure 
This is a modification of 3.5.2. where a parameter a 
varying as [0,1] is considered 


| vYty)dy 
Ya) 


guie (18) 


Y (y) dy 
s Yta) 

The idea which underlies this process is that the 
threshold level a permits an estimation of yq taking into 
account that elements which have grades of membership 
lower than a do not need to be considered. In other 
words, values (regarded as) non-significant can be 
eliminated from the estimating process. 


3.6  Alteration/Correction of the 


calculated values 

Once the defuzzification stage is completed the 
controller's implementation may also be regarded as 
achieved. Two alternatives can, thenceforth, frame the 
designer's work: a) or the results already obtained 
satisfy his/her endeavoured objectives and then no 
further work is needed, or b) on the contrary (because 
unsatisfactory), such results require a reappraisal of 
some of the steps previously performed. 

In the latter case, following an approaching strategy 
that proceeds from the easier to the more difficult 
alterations which the whole design may undergo, the 
designer has at his/her disposal several other 
possiblities which may be briefly sketched as follows. 

A) Introduction of slight variations in the parameters 
of the curves which initially described the selected 
linguistic attributes. 

B) Alteration of these attributes and/or introduction of 
new ones. 

C) Alteration either of the range of the intervals into 
which the universes ol discourse had been divided or 
of proper universes. 

D) Introduction ofone (at least) new input variable, 

Typically, if the designer has, at the outset, 
considered one input variable only (usually 
corresponding to that "error" which is given by the 
difference between the desired state and the actual 
values provided by the controller) then a good 
strategy consists in the further inclusion of an "error 
variation" working as a second input. Obviously, 
this will lead to the substitution of (2) by more 
complex inference rules given by (3), to the 
specification of its universe of discourse, to the 
definition of new linguistic values and new 
membership functions for Lhe error variations, etc. 

E) Substitution of the classical operations (union, 
intersection, direct product, max-min composition, 
etc.) by more elaborate operations such as those 
exemplified in the t-norms or t-conorms previously 
referred to. Different values of the parameters w or À 
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there focused must, in these conditions, be also 
taken into account. 

HF) Construction of a more elaborate controller where 
"adaptation" and "learning" have, from the outset, 
been already included. 

Insofar as we know step É) has not yet been the 
objective of any specific research. On the contrary, step 
F) is already used in many of the present-day controllers 
but its underlying theory is outside the scope of this 
work, About steps A) to D), the most frequently used in 
the literature of fuzzy controllers, the case-study we 
present in the sequel is clear enough to provide the 
reader with a description of the whole design philosophy 
previously referred to. 


4 An Illustrative Example 

This case-study treats a real, concrete problem which, 
presently, is the objective of intensive research: the 
construction of a fuzzy engine which can control a third 
order, non-linear system with ill-defined, time-variable, 
delays. This system consists, broadly speaking, of an 
electric motor, a gear-box and a load-torque disturbance 
(represented by an iron-bar). The objective of the control 
is that the motor holds the bar in some pre-defined 
equilibrium position. 

The whole system may be graphically displayed as in 
fig.7. 


LDAD-TORQUE 


DISTURBANCE 


att 


HOTOR 
GEARS 
LOAD ] 


fig. 7 


At the outset, the fuzzy controller was regarded as a 
SISO system (albeit subjected to the load disturbance 
Tp) where Lhe input variable is an 'error' E=9(t) - O 
and the output (representing the control variable) is the 
tension V applied to the motor. 

In the following considerations, and for the sake of 
simplicity, the universes of discourse of these variables 
will already be (assumed as) normalized. Thus, if the 
value O is supposed to be the representative of the 
desired equilibrium position, then to the limiting 
boundaries of such universes will correspond normalized 
intervals |-1, +11. 

Two types of fuzzy predicates were, initially, assumed 
to describe Lhe possible variations of E and V. "Small" 
(S), "Medium" (M) and "Great" (G) for E and "Decrease" 
(D), "Null" (N) and "Increase" (1) for V. Their 
membership functions obey the representation a) of 
fig.4. 

At the outset, the relationships between these 
predicates (nothing but the interference rules 
implemented in the fuzzy engine) were supposed to obey 
the fuzzy propositions 


R,: “IF E=M THEN V=Nº 
Ro:"IF E=G THEN V=Dº 
R3:"IF E=S THEN V=P 


Thus, assuming for the sake of exemplification that 
the universes of discourse of É and V were discretized in 
five Intervals (in reality such yniverses comprise 1000 
intervals), the partial matrix which expresses Rj would 
be given by means of the inference operation "x" (which, 
in this case, corresponds to the t-norm b) of Section 3.3.1) 


11x [0002081]=[000208 1 
0.8 000208 08 
0.1 0001011 
0 000 0 01 
0 000 00 


and the global input/output relationship Ry (obtained 
using the usual union) would be numerically described 
by 


O O 0.2 0.8 1 
O 0.2 0.2 0.8 0.8 


Ro =/0.1 0.3 1 0.3 0.1] 


0.8 .0.8 0.2 0.2 O 


| 0802 0 O 
As a further illustrative point, if we assume that 
8=0.4 then to this value would correspond a concrete 
defuzzified value V=0.52. The operations which underly 
the evaluation of this value are graphically displayed in 
fig.8, quite similar to those which fig .6 expresses. 


V = "Null 


E = "Great” V = "Decrease” 
fig. B 

As a final, albeit provisional result, the overall 
response of the system is that depicted in fig.9. 

It can be observed that this response undergoes some 
severe disadvantages, namely the existence of an 
undesirable overshoot and, above all, an undesirable 
steady error (for a disturbance which, in the real case, 
corresponds to a step whose amplitude is 8 N/m). 
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fig. 9 


In order to overcome these problems a new design 
strategy was implemented. The new approach has laid 
emphasis on two main types of alterations: the inclusion 
of an aditional input variable (the “error variation' 
DE=68 (t+1)- O (t)) and the introduction of new fuzzy 
predicates for E and V. 

From this prespective E is now described by means of 
seven fuzzy predicates ("positive”, "small positive”, “zero 
positive", "zero", “zero negative”, "small negative”, 
“negative”) while the predicates for the output V become 
“increases very very much", "increases very much”, 
increases", "does not act”, "decreases”, "decreases very 
much", "decreases very very much", About DE, its 
predicates correspond to "great positive”, “positive”, 
“null”, "negative", "preat negative”. 

For the sake of simplicity we shall not present neither 
the 35 rules which the combination of the predicates 
leads to nor the calculi they entail. The unique aspect 
which the reader must, however, bear in mind is that 


TFRE=. AND DE=.o THEN V= o" 


where the "AND" was assumed to correspond to the 
usual fuzzy intersection. 

As a consequence of these alterations the response of 
the system is now that one depicted in fig. 10 where the 
inconvenients that the former situation encompassed 
have already been overcome. 


! 


DISTURBANCE 


5 Outline 

Linguistic control is a recent and unconventional 
control technique which overcomes a great deal of the 
traditional problems that classical control is still facing 
nowadays. In this paper the main stages on which the 
design of a general linguistic controller is based were 
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briefly examined. No reference was, however, 
purposively made to the analysis of the statical and 
dynamical properties of the controller (beseting several 
indexes which deal, for example, with the completeness 
of the rules, their possible interaction, their consistency, 
etc.). Likewise, the robustness and stability of the 
overall controller are matters which, due to their 
complexity, require so detailed an approach that they 
are clearly out of the scope of the present analysis. 
Despite these restrictions (the examination of which 
will, none the less, be deferred to further works) the 
authors believe that the foregoing exposition is able not 
only of providing a self-content theoretical and pratical 
tool on which the understanding of a great many deal of 
papers on fuzzy control becomes thenceforth possible but 
(more important, perhaps) also of stimulating the reader 
for the exploitation of a subject which is full of research 
promises, 

As Kaufmann [6] asserts: "...At the present we are 

unable to design machines that can compete with 
humans in the performance of such tasks as recognition 
of speech, translation of languages, comprehension of 
meaning, abstraction and pgeneralization, 
decision-making under uncertainty and, above all, 
summarization of information. In large measure, our 
inability to design such machines stems from a 
fundamental difference between human intelligence, on 
the one hand, and machine intelligence, on the other. The 
difference in question lies in the ability of the human 
brain - an ability which present-day digital computers do 
not possess - to think and reason in imprecise, 
non-quantitative, fuzzy terms. Itis this ability that makes 
possible to decipher sloppy hand-writting, understand 
distorted speech and focus on that information which is 
relevant to a decision. Anditis the lack of this ability that 
makes even the most sophisticated large scale computers 
incapable of communicating with humans in natural - 
rather than artificially - constructed languages. 
“In our quest for precision we have attempted to fit the 
real world to mathematical models that make no 
provision for fuzziness. We have tried to describe the laws 
governing the behaviour of humans, both singly and in 
groups, in mathematical terms similar to those employed 
in the analysis of inanimate systems. This, in my view, 
has been and continue to be, a misdirected effort, 
comparable to our long-forgotten search for the 
perpetuum mobile and the philosopher's stone. What we 
need is a new point of view, a new body of concepts and 
techniques in which fuzziness is accepted as an 
all.pervasive reality of human existence. Such concepts 
and techniques could open many new frontiers in 
psychology, sociology, political science, philosophy, 
physiology, economics, linguistics, operations research, 
management science and other fields, and provide a basis 
for the design of systems far superior in artificial 
intelligence to those we can conceive today." 

That such hopes soon become realities is, the last but 
not the least, what the authors expect from the readers 
ofthis work. 
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RESUMO 


No contexto dos processos digitais cicloestacionários 
abordados neste artigo, mostra-se que as cadeias de 
Markov podem ser um instrumento valioso para a 
modelização dos sinais e o cálculo de espectros de 
potência. 


1 Introdução 

Em muitos sistemas de comunicação digital, a 
informação é transmitida sob a forma de processos 
aleatórios síncronos. À estacionariedade das sequências 
de símbolos incorporadas nestes processos - implicando a 
variação periódica, no tempo, dos respectivos 
parâmetros estatísticos - permite defini-los como 
processos cicloestacionários [1]. A aplicação do conceito 
de cicloestacionariedade é vantajosa para fins de análise 
espectral e pode desempenhar papel importante na 
concepção e avaliação de circuitos de sincronismo [2-4]. 

Com frequência, os processos cicloestacionários 
usados em comunicação digital podem ser descritos por 
modelos markovianos. Estes modelos são 
particularmente úteis na caracterização de modulações 
“com memória” [5]; uma adequada modelização deste 
tipo permite, nomeadamente, aliar rigor e eficiência 
computacional no cálculo de densidades espectrais de 
potência [5-8]. 

No presente artigo, abordamos a aplicação de um 
modelo markoviano ao estudo de processos digitais 
cicloestacionários: uma caracterização genérica destes 
processos é o objectivo da Sec. 2; à aplicabilidade do 
modelo markoviano é discutida na Sec. 3, 
exemplificando-se na Subsec. 3.2 para um formato de 
modulação "com memória" de grande generalidade e 
interessse prático. 


2 Processos Digitais Cicloestacionários 


2.1 Caracterização genérica 

Iremos considerar neste artigo sinais aleatórios 
associados a sequências de símbolos pertencentes a um 
alfabeto finito e redutíveis à forma: 


+ uu 
s(t) = > zU-nT). (1) 


n= =» 


Original submetido a publicação em 17 de Setembro de 199] 


ABSTRACT 


This paper deals with ciclostationary digital 
processes. In this context, a Markov chain approach is 
shown to provide a powerful tool for signal modelling 
and calculation of power spectra. 


p,= Prob (e, (8) = w (0) (2a) 


(2b) 


Pie = Probte  (=w (Oz (D=w(0)) (mz1) 

Todos os possíveis zn(t) são sinais de energia 
(admitiremos, ao longo deste artigo, que são funções 
reais de t), e, para cada n, o “impulso” zyíLt) é 
seleccionado, no âmbito de um mesmo conjunto de M 
impulsos possíveis (wi(t), wo(t),...., wm(t); (M finito), de 
acordo com a sequência de simbolos e uma dada regra de 
codificação. Admite-se sempre que a sequência de 
impulsos exibe estacionariedade (no sentido lato), não 
dependendo de n, por conseguinte, as probabilidades 
seguintes: 

e 
m (t) = E lb] = > Elz U-nT 9] 


n=-—-m 


R, (t+)=E Is(t) stt+ )]= 
+ DD 


— Ny ) — , — | 
>. 4 Elz, (t n; Pa Rio no dl, 


n,=-0n,=-D I 


] 2 
É fácil verificar que, nestas condições, os processos 
são cicloestacionários (no sentido lato), pois a média 
ms(t) e a autocorrelação Rs(t, t+), dadas porsão funções 
periódicas de t (período T,): qualquer que seja j inteiro, 
ma(t+]f;) = me RL +jTs, t+H)T;+u) = Ret, L+v). 
Invocando (2a) e (3a), temos o sinal periódico 


& 


J 


m dB = >: w (t-nT) (4) 
em que 
.  M 
Ww = S P, w LL) (5) 


Este sinal periódico pode ser escrito na forma 
+ ma 

m(y= > 
É) = 


n=-w 


c exp Q2nnuT ), 


sendo 
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TR + mo 
l 5 = us 
c =—| / wtt=n7 9 expt=;2nnt jdt= 
à TI=PBb... s 5 
5 8 = — 0 
| + o. — | 
= — | w(t) expt—y2nntT. PE w(— 5 (6) 
Eq — DO T T 
8 g E 
com Wí(fj= Flw(t)] (F: transformada de Fourier; 


C.n=Cn*). Naturalmente que, atendendo a (5), 
M 
Wp=>DpWi), cmWiA=Flw (o). 
1=1 
O sinal s(t)=s(t)-m,(t) pode ser escrito na forma 
ju 
s0= 5 AM =87.) (7) 
C de a 
n=-uw 
com impulsos z,'cMt) seleccionados no âmbito do conjunto 
IwylcXL), woteMt), ..., wmicHL)y, em que wjteML) = wi(L)-w(L), 


jJ=1,2,..., M. Obviamente que são iguais a pj e p;k!M) 
(ver egs (2)), respectivamente, as probabilidades 
| (8a) 
Prob tz 0) =w “(t)] 
(8b) 


te) o teste. Ve) | 
Prob Es mb=W, (uz, ()=uw, (0) (mz1) 


Note-se que stcMt), tal como s(t), é cicloestacionário de 
período T,, e E [steM(t)]J=0. 


2.2 Cálculo de espectros de potência 

Sendo s(t)=mu(t) + st), a densidade espectral de 
potência de s(t), Gu(f), apresenta duas parcelas: uma 
parcela de espectro discreto (espectro "de riscas"), 
Gs!df), correspondente à componente periódica mt), e 
outra de espectro contínuo, Gçtc)(f), correspondente à 
componente st). 

Atendendo a (6), conclui-se facilmente que 


di À ie Sds É (9) 
Gb= SIW > DT; =. 


= — Bb 
s Fl 5 


Para a determinação do espectro contínuo há que, 
atendendo à cicloestacionariedade do processo, começar 
por achar a média em t de E I|sdt),s(t+1)], num 
qualquer intervalo de duração T,. A transformada de 
Fourier da função de autocor rálação obtida (verl9]) é a 
densidade espectral de potência 

+ qu 


dimmo 


BRT) exp (-y2n/mT ) (10) 


E, m=-u 
em que Em() = ElZn + m! HDL te D) CD) = Flzj(1)]). A 
estacionariedade admitida para a sequência de impulsos 
implica que E .n(f) = Em *(f), pelo que 
Pn=E o +p de] É e E mt expt=p2nfmT | (11) 
é E = 
O cálculo de Guto(f) é fácil recorrendo às 
probabilidades p; e p; k(m), que permitem obter 
S (129) 
ni E tel a 
EqD = Le Iwo] 


t=1 


0-5 E nte Wo Wen mer (12 


i=1k=1 


(ein = ogia] efa 50) 


Casos Pa de EE 


a) Na ausência de correlação entre os sucessivos 
impulsos zn(t-nT.) (modulação "sem memória”), 
podemos dizer que p; k'M)=pk, necessariamente, 
Nesse caso temos, para mz, 


se 
e wie) 
Et = o ep >» pre 

is 
e portanto o espectro anil reduz-se à primeira 
parcela de (11), com Eo(f) dado por (12a) (Aliás, se 
houver correlação mas apenas até um dado 
Mmáx 20, 


m 
+ 00 max 


torna-se possivel substituir > por Ea em (11), 
m=1 m=1 
por ser nesse caso E (f)=0 para m> Máx). 
Por outro lado, o espectro contínuo pode também 
ser apresentado na forma seguinte: 


M 
l 
G'Hp=— Sos p IwS Hp- win = 
8 T , k 
gs 1=1 k>1 
M M 
| SS” | 
E s p > plWin-W (pI (13) 
op SR | k 
s |= k>1 


(ver em [10], a propósito, o cálculo da densidade 
espectral de potência para PPM digital). 

b) Nas situações de "modulação linear”, isto é, quando 
Zn(t)=anw(L) e wi(t)=abw(t), j=1,2,..., M (ou seja: 
para cada n, a, é seleccionado no “Ambito do 
conjunto (a!!), at2),..., at(M))), o espectro discreto é 


+ q 
Pivaê > do) 04 


l 
p= 2 ( 


M 
coma = Ela |= > pao. 
É Bê 
O espectro contínuo pode ser obtido a partir das egs 
(11)e(12), SME 


(15a) 
E) 2 IWipp 


(15b) 


E (A IwWpP,mz1 


r R im) (a? 


M 
= Ps a (DD or 
=Ela ]-(a) = > pia | —ta) 


1=1 


com della -— 1º] 
ã n 


M M 
| > a oi Um) (4) 4h) 
e, para mz1,R tm) Ela a im! 22 Pi Fr US. 


Repare-se, combinando as eisstica (10), (14) e (15), 


que o espectro de potência pode ser escrito na forma bem 
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conhecida [11] 


+ o 
> R (m) exp ( =J2nfmT ) (16) 


É ad 
G = q IWA) | 

Emi 
em que R,(-m)= Elagan.mlI= Ralm) e Rj(0)= Elaç2]= 02, 
+ (a)2, isto é, o espectro de potência do sinal PAM 
assente numa dada sequência e num dado impulso 
elementar é produto de dois factores, um dos quais é o 
espectro de potência da sequência e o outro depende 
unicamente do impulso elementar (mais precisamente, 
da sua densidade espectral de energia). 

Note-se ainda que, no caso linear, a ausência de 
correlação significaria que, para mz 1, Ry(m)=(a)? - ou 
seja, Em()=0- e portanto 

: Eai 
G (= T o WO) (17) 


3 Aplicação de um Modelo Markoviano 


3.1 Características gerais do modelo 

O emprego de cadeias de Markov na modelização de 
processos digitais ciclvestacionários tem sido proposto 
por diversos autores, designadamente nas referências [5- 
8). As cadeias de Markov usadas nesta modelização 
referem-se concretamente a sequências de “estados” 
associadas aos sinais, e são, de um modo geral, cadeias 
homogéneas, isto é, as probabilidades das transições de 
estado para estado são estacionárias [5]. Sendo T, o 
período dos processos cicloestacionários, é conveniente 
definir um modelo markoviano em que estas transições 
se processam de T, em T,; a cada sequência de formas de 
onda como a genericamente indicada na relação (1) 
corresponde assim uma sequência de estados de ritmo 
idêntico, ambas de acordo com a sequência de símbolos 
incorporada no sinal e com uma dada regra de 
codificação. Numa modulação "com memória”, em que o 
emprego do modelo markoviano é vantajoso, o estado 0, 
atingido no instante de ordem n depende, não só do 
símbolo ou do grupo de símbolos “correntes”, mas 
também do estado atingido no instante precedente, de 
ordem n-1: a dependência em relação à história passada 
do processo faz-se unicamente através de 0,1. Às cadeias 
de Markov são assim completamente caracterizadas 
pelas transições, ditas de "primeiro grau”, entre o estado 
num instante e o estado no instante seguinte. Numa 
cadeia de Markov homogénea com L estados, as 
probabilidades das transições de grau m 

a = Prob (0, mto, = Uimel;Lk=1,2... o 
podem, para qualquer m, ser obtidas unicamente em 
função das probabilidades P, !'!!, no seu conjunto 
apresentáveis na forma de uma “matriz de transições" P 
com L linhas e L colunas. 

Admite-se que, nas aplicações correntes, as cadeias de 
Markov são "completamente regulares" [5] e assim: 


a) o conjunto das probabilidades estacionárias dos 
diversos estados, Pi (i=1,2,..., L), constitui solução 
única da equação matricial [Pl P2 ..PL|xP=[P1 
P2....PL) 
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L 
(seimpusermos > P=1eP 20,i=1,2,..,L). 
i=1 
b) Verifica-se que Pi,k(m) é o elemento da linha i e da 
coluna k da matriz Pm. 


c) 
PE, 
PP; 
Em PPS asma (19) 
m—+ 
PPP, 


Casos há em que se verifica uma correspondência 
biunívoca entre estados e formas de onda possíveis |[6,7|, 
isto é, L=Me 
2 t=uwleo = 20) 

Naturalmente que, nestes casos, e atendendo às 
relações (2) e (20), se tem 
p= Pp (21a) 
(im) — ptm) (21b) 

Pik Ri Pa 

Em [5] e |8] são abordadas situações em que o modelo 
markoviano não pressupõe correspondência biunivoca 
entre os estados e as formas de onda. O chamado "código 
bipolar" é um caso concreto em que é possível definir um 
modelo markoviano em conformidade com os princípios 
enunciados mas sem correspondência estados / formas de 
onda (ver Apêndice A). 


3.2 Modelo markoviano com 
funções-estado 

O modelo markoviano com funções-estado é aplicável 
em muitas situações correntes (ver [7] para um 
tratamento detalhado e bastante geral). Nestes casos 
qualquer dos ramos do diagrama de transições de estado 


é do tipo representado na figura 1, pelo que ná 


ki 


figura 1 


correspondência estados/formas de onda e se verificam 
as relações (21). Isto permite, nomeadamente, um 
cálculo sistemático das densidades espectrais de 
potência que concretiza, com simplicidade, o exposto na 
Sec.2: uma vez obtidas, e organizadas numa matriz p, as 
probabilidades p; k'l), podem as probabilidades p;k'"! 
(m>1) ser obtidas directamente através do cálculo de 
pm 

Em |12] é apresentada uma aplicação do modelo 
markoviano com funções-estado a sinais “de tipo 
OQPSK" (de tipo MSK) e outros relacionados com estes 
(mais precisamente, uma aplicação às respectivas 
componentes em banda de base). Tira-se aí partido dum 
formato generalizado para a representação de sinais 
daquele tipo (ver, a propósito, [13,14] e o Apêndice B): 
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descrito o sinal de tipo OGQPSK sob a forma de um 
esquema de codificação ENCAP-4N, o número de 
funções-estado é de 22N+3 ocorrendo a ritmo 1/T,, com 
T.=2T: (T: tempo de duração de um bit). A precisão do 
método de cálculo de espectros por esta via é, 
evidentemente, elevada, desde que a precisão do cálculo 
das 22N+3 transformadas de Fourier Wi(f) também o 
seja. É possível, todavia, aumentar a eficiência 
computacional sem afectar a precisão do cálculo, 
bastando explorar o facto de diferentes sequências de 
bits corresponderem a funções idênticas ou simétricas. 
Em [13] mostra-se que a densidade espectral de potência 
das componentes de um qualquer sinal ENCAP-4N pode 
ser calculada na base de um modelo markoviano com 
22N+2 funções-estado, e como são simétricas duas a 
duas, há apenas que calcular Ng=22N+1 funções Wi(f) 
(admite-se a transmissão duma sequência (vn) de 
simbolos binários estatisticamente independentes, sendo 
vn=+1 ou -1 com igual probabilidade 1/2). Chega-se à 
fórmula 


“o | "o "o R 
Gip=dAp= — | Ny WD" +2Re| > + WOW. 
: Ú 1=1 p=1k=1 (22) 
N 
> (2% — PN soa) espt-sanfmro || 
m=1 U 


As funções wilt) (1= 1,2,... Ny) correspondem, para N 
impar, aos diversos casos particulares de LI, q(L)+ Is (Lt), 
com v.n.j=1. O indice é 


AN 
i=14D,4 5 2rc (23) 
1) a. nu 
n=i 
Para N par, as funções wilt) (i=l, 2, No) 


correspondem aos diversos casos particulares possíveis 
de Qu o(t)+ Qua(t), com vn=1. O índice é igualmente 
dado por (23). 

Quanto ao cálculo das probabilidades p, pm), a 
principal dificuldade é para m=1, como se depreende do 
que vimos atrás, pois há que ter em conta a estrutura do 
sinal, traduzida nas regras de codificação do Apêndice B 
(que um diagrama de transições de estado poderia 
ilustrar). Para m>1, basta calcular sucessivas potências 
da matriz das transições, Às probabilidades p, '!! podem 
ser obtidas aplicando as seguintes regras: 

a) Seja N=2L+1, com L=z0. Quando à for dado por 

(23), se CN 9=0, então 
l/4, para k= | + 4D, - 


2N - 2 
(1) ", n+Z | 
= A = F | | 
k + DLC j=1238 (24) 
n=1 
O, para os restantes valores de k, 
E se Cno = |, então 
l/4, para k=2N + | =.) 4D, 6 — 
2N —2 
(1) s, n+2.. 
| E E , ++ ] - 2: 
Pik dá dia nO e do (25) 
n=1 
O, para os restantes valvres de k, 
26 


b) Seja N=2M, com Mz1. Quando i for dado por (23), 
se Cn.19=0, então p; k!'!' é dado por (24), e, se Cn. 
10=1, então p; k'! é dado por (25). 

c) Tendo calculado metade dos elementos da matriz 
das transições, acha-se facilmente a outra metade, 
com base na relação 


1) SA qtos 

Pont 1-1,2NG+1-A Ph y Lk=1,2,...,2N, (26) 
que aliás, arrasta esta outra, para m> 1: 

tmi AME, à Er é 
Pongti-12Nç+I-k Pi Wk= Loo (27) 


(consultar [13] para detalhes). 


Note-se, a terminar, que a modelização usada em [12] 
tem a vantagem de poder ser usada tanto pura as 
componentes "em fase" e "em quadratura” de um sinal 
de tipo OQPSK, s(t), como para as componentes dos 
sinais st), s&1), etc (cadeias de Markov similares, mas 
com diferentes conjuntos de 22N+3 funções-estados), o 
que se reveste de bastante interesse na análise dos 
subsistemas de sincronismo [12,15]. 


4 Observações Finais 

No presente artigo, mostrou-se que a aplicação de 
cadeias de Markov na modelização de processos digitais 
cicloestacionários pode ser bastante vantajosa, U modelo 
markoviano permite, nomeadamente, sistematizar a 
determinação das probabilidades p; k'M', associadas às 
transições entre formas de onda, o que se reveste de 
muita utilidade para o cálculo rigoroso dos espectros de 
potência. À bem conhecida fórmula de [6], para o cálculo 
desses espectros em casos em que o modelo markoviano 
joga com funções-estado (muitas vezes mal citada, talvez 
devido a erros Lipográficos presentes nesse relatório), é 
justificada neste artigo com relativa simplicidade, no 
quadro de uma abordagem bastante geral que, 
consentindo a aplicação do modelo markoviano, não o 
pressupõe necessariamente. Neste sentido, a abordagem 
feita em [7] para esses casos difere da aqui adoptada, 
algo mais simples e explorando directamente a hipótese 
de cicloestacionariedade dos processos. 

Como exemplo de aplicação das cadeias de Markov, é 
indicado um método prático para obtenção dos espectros 
de potência de uma vasta classe de sinais “de tipo 
OQPSK” [14], incluindo muitos de utilização corrente 
(modulações MSK, GMSK, TFM, IJF-OQPSK, etc). 


Apêndice A 
Código bipolar 


Suporemos aqui que os simbolos da sequência da 
entrada são independentes e que Prob(1)= Prob(0)= 1/2. 
Aos simbolos 1 correspondem alternadamente os 
impulsos wi(L) e wo(Lt), e aos símbolos O correspondem 
sempre impulsos wa(t), com 


vw (0) = w(t) (A.1a) 
wo) = — wlt) (A 1b) 
ww LL) = 0 (A 1c) 


podendo o impulso elementar wít) ser o impulso 
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rectangular padrão, com duração T,, mas não 
necessariamente. Para este código de linha, a memória 
introduzida no processo traduz-se num modelo 
markoviano (cadeia homogénea e completamente 
regular) com L=2 estados e o diagrama de transições 
representado na Fig A.l. 


1/wy 


1/w, 


figura A.1 


Atendendo à figura, e de acordo com o exposto na Sec, 
3.1, verifica-se que, neste caso, a matriz de transições 
tem os quatro elementos todos iguais a 1/2 (note-se que 
Py !D=Po 9! =Prob(0) e P; g!=Po !b=Prob(1)), pelo 
que, também para m>1, Pjfmi=1/2 (1]= 1 ou 2),e 
ainda Pj=P,=1/2. 

Atendendo à figura A.1, as probabilidades p; tm) 
(1,k=1, 2 ou 3) das transições entre as formas de onda 
wi(t), wo(t) e wa(L) podem agrupar-se em matrizes que de 
alguma forma se relacionam com as matrizes Pm (Por 
exemplo: pra! =0,p21! Li= P, 9! L), py atb= PP; 9!D), etc.: 
e, para m>0, pratm+ti)=Po tmP, otb, 
po (m+ 1)= P, tm)Pj 3! | , pa tm+ 1) = PjPy tmiP, a! 1) 
+ PoPo jtm)P; ot) etc.). Então, 


O 1/2 1/2 
pal = 12 0 1/2 (4.2) 
1/4 1/4 1/2 
e, param >0, 
1/4 1/4 1/2 
Lp; = 14 1/4 12 (A.3) 
1/4 1/4 1/2 


sendo ainda notório que p;=1/4, po=1/4 e py=1/2. 
(Note-se que [p;u'm'|zl|p;klm: a cadeia de Markov 
refere-se aos dois estados (1 e 2), e não às três formas de 
onda). 

Obtidas as probabilidades p; e p; k'M), torna-se fácil o 
cálculo do espectro de potência com o método geral 
descrito na Sec.2, particularizado para a situação de 


linearidade. Neste caso, em que alDl=1, al2)=-1 e 
atyi= (0), temos 
— A 4 
a=0 A 
R (0)=0/=1/2 Voa 


A. Gusmão 


TA ie A 


Em) di) 4Ã) 
E aa ads: 


O ,m>1 (AM) 

O facto de termos a =0 implica a ausência de espectro 
discreto, sendo o espectro continuo obtido combinando as 
relações (16) e (A.4): 


l 
E RA PS ' 
G, (N=G (= 2T [1-5 Qnfro |) (A 5) 
8 


Apêndice B 


Os sinais do tipo OPSK podem ser descritos na forma 


s()= | [o +IQ do) | explj2nf + 05)| 


em que 


m 
= yr E 
IO=D1 dd) e QU=20,[8 
Fm 
sendo as funções Is mtt) e Qsm(t) nulas fora do intervalo 
[mT,(m+ 1)T]. De acordo com as regras de codificação 
ENCAP-A4N [13]: 


= UNO pjtt— 2n7), 


Q, ant) El nYplt2n +9)T-+]; 


E on = o qt DF bi 
O, on 10=U9,- nYmlt— (+ 17], 


Os índices das funções y são obtidos, com base na 
sequência (vn), da seguinte forma: 


2N À 
fo » gP-lo (B.la) 
SS da fi 
P= 
2N 
k= 2 2", (B.1b) 
p=1 
2N 
l= E 92N — pp (B.1c) 
pn 
p=1 
2N 
-= “º oP-l 
m=2 2""D, (B.1d) 
p=1 
com 
É (B.2a) 
Ca Van anna p=12,...2N) 
Dm hr Cgi | 68.26) 


(Vj=1 ou O consoante v;= +1 ou -1; & significa “adição 
módulo 2") 
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RESUMO 


Quando os sistemas passivos de condicionamento 
ambiental se revelam insuficientes na obtenção do 
conforto térmico dos edifícios, ou quando o nível e tipo 
de consumos de energia nos edifícios exigem o recurso 
a equipamentos energéticos cuja eficiência de con- 
versão da energia primária é normalmente próxima 
dos 35%, a bomba de calor poderá constituir, em certas 
condições, uma tecnologia alternativa no uso da 
energia. Com efeito a bomba de calor é um sistema que 
permite aumentar a eficiência do uso da energia 
sempre que se pretende dispôr de calor a baixa ou 
média temperatura e sempre que a sua eficiência 
(COP) exceda o valor 3. A sua economicidade é 
particularmente relevante nas situações em que se 
requerem refrigeração e aquecimento simultâneos, 
Estes sistemas têm pois amplo âmbito de aplicação nos 
edifícios, muito particularmente nos sectores não 
domésticos, onde é grande o nível e a diversidade de 
consumos de energia pelos vários tipos de utilização. 

Com este trabalho pretende-se ilustrar a aplicação 
da bomba de calor aos edifícios no que respeita a 
desempenhar com vantagem algumas funções, 
designadamente: recuperação do calor degradado, 
produção combinada de frio e calor e recuperação do 
calor evacuado de uma piscina, 


l Introdução 

Os consumos de energia nos edifícios representam 
entre nós apenas cerca de 20% do consumo energético 
nacional, devido à amenidade do clima por um lado e ao 
baixo nível de conforto térmico por outro, que 
genericamente e tradicionalmente caracteriza os 
edifícios em Portugal. O aumento dos custos de energia 
primária, as crises energéticas, a existência de uma forte 
taxa de crescimento dos consumos nos edifícios, 
resultante da crescente exigência de conforto por parte 
da população, de deficientes projectos dos edifícios do 
ponto de vista energético, de sobredimensionamento e/ou 
inadequação dos equipamentos e da deficiente regulação 
e controlo dos sistemas, conduziram a uma nova filosofia 
que é vulgarmente designada por Utilização Racional de 
Energia. 

O condicionamento passivo dos edifícios e a utilização 
da bomba de calor constituem exemplos de alguns 
métodos de utilização da energia que se enquadram 
nesta filosofia. Este último exemplo poderá ser encarado 
quando haja situações que não permitam dispensar 
sistemas energéticos para conferir aos espaços as 
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ABSTRACT 


The heat pump becomes an alternative technology 
in the use of energy in buildings whenever passive 
technologies of environmental conditioning are not 
enough to satisfy thermal comfort or the energy 
demand requires the use of conventional equipment 
converting primary energy with an efficiency of only 
35%. Heat pumps offer considerable scope for 
improving the efficiency of energy use wherever heat 
is required at low temperature (below 100ºC) and the 
thermodynamic efficiency (COP) of the heat pump 
exceeds the value 3. 

Their application is viable in buildings, 
particularly in the non domestic sector, where a 
variety of simultaneous energy uses is needed. Its 
economic advantage can become particularly 
important in the combined heating and cooling 
functions. 

The analysis of some heat pump functions in 
building is developed in this paper, namely the 
recovery of rejected heat, the combined cooling and 
heating and the combined heating and 
dehumidification effects. 


condições de conforto ou quando haja outros tipos de 
utilização em que se pretenda dispôr do calor a baixa ou 
média temperatura. Neste trabalho, pretende-se 
explorar as potencialidades da bomba de calor no que se 
refere à sua aplicação aos edifícios. São analisados 
alguns exemplos de funções que a bomba de calor pode 
desempenhar, com vantagem, num edifício. 


2 Consumos de Energia num Edifício 
por Tipo de Utilização 


2.1 Sector doméstico 

Como se distribuem nos edifícios os consumos de 
energia pelos vários tipos de utilização”? 

Segundo um trabalho recentemente publicado pela 
Direcção Geral de Energia a distribuição percentual dos 
consumos de energia final' no sector doméstico pelos 
vários tipos de utilização, correspondente ao periodo de 
12 meses anteriores a Junho de 1988, foram: 


1 Energia final - energia disponivel no local de utilização sob a 
forma utilizável. Ex: electricidade nas tomadas, gasóleo nas 
caldeiras, gás na cozinha, água quente no termo-acumulador, etc. 
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Água quente sanitária e cozinha 78% 
Electricidade específica? 12% 
Aquecimento ambiente 10% 


As formas de energia mais consumidas para o 
aquecimento ambiente são a electricidade e a lenha e as 
mais utilizadas para a água quente sanitária e cozinha 
são o gás em garrafa e /ou canalizado e a lenha. Todavia, 
é razoável admitir, que face à crescente exigência de 
conforto nos edifícios, cresça o nível e se altere a curto 
prazo esta distribuição de consumos e que até surjam 
novos tipos de consumo (p. ex. arrefecimento ambiental). 

Uma estimativa das necessidades de energia para 
uma habitação do tipo T4, com inércia moderada, 
ocupada por agregado familiar de 4 pessoas, aponta para 
a seguinte distribuição percentual de consumos médios 
anuais de energia útil 3. 


Aquecimento ambiente 37% 
Arrefecimento ambiente 31% 
Água quente sanitária e de cozinha 23% 
Electricidade específica 9% 


Nestes cálculos consideraram-se os seguintes 
parâmetros de referência: 


Aquecimento ambiente 


Temperatura interior do projecto 20“C 
Temperatura exterior do projecto 6º 
Duração do período de utilização 900 h/ano 
Potência correspondente à ponta 
de carga térmica 70 W/m? 
Arrefecimento ambiente 
Temperatura interior do projecto 24º€ 
Temperatura exterior do projecto 32ºC 
Duração do período de utilização 540 h/ano 
Potência correspondente à ponta 
de carga térmica 100 W/m2 
Água quente sanitária e de cozinha 
Consumo de água 200 I/dia 
Temperatura da água quente 45"C 
Temperatura da água da rede 10ºC 
Duração do período de utilização 365 dias/ano 


2.2 Sectores residenciais não 


domésticos 


Para estes sectores a distribuição dos consumos é 
muito variável devido à sua heterogeneidade no que se 
refere às funções e serviços a satisfazer. 

Neles se incluem os hóteis, escolas, escritórios, 
hospitais, pavilhões gimnodesportivos, entre outros. 
Todavia será razoável admitir que o consumo de energia 
para a climatização ambiental (aquecimento/ arrefeci- 
mento) deverá ser o que mais peso relativo tem no 
consumo total de energia nestes sectores. 


2 Electricidade específica - energia eléctrica consumida em 
equipamentos especificamente eléctricos (Ex: ferro de engomar, 
aspirador, máquina de lavar, frigorífico, etc.) para os quais já não 
são utilizados equipamentos alternativos consumindo outra 
forma de energia. 

3 Energia útil - energia efectivamente utilizada para salisfazer 
uma dada necessidade e que resulta do processo de conversão 
final. Ex: água quente utilizada, calor consumido no fogão, calor 
para conforto ambiente, etc. 
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Deverão ainda figurar nestes sectores novos tipos de 
consumo que não têm expressão no sector doméstico, tais 
como: 


- Renovação do ar em edifícios com taxas de ocupação 
muito elevada ou com exigências específicas de 
qualidade do ar; 

- Desumidificação de ambiente com taxas elevadas de 
produção de vapor; 

- Produção combinada de frio e calor. 


Face à grande diversidade de consumos e à exigência 
cada vez maior de conforto nos edificios sem dispêndio 
excessivo de energia, que soluções técnicas alternativas 
se oferecem? 


- Recurso a fontes energéticas gratuitas: solar, calor 
ambiente, calor residual de eventuais efluentes 
térmicos, 

- Melhoria dos equipamentos energéticos 
convencionais (caldeiras, ar condicionado, painéis 
radiantes, permutadores de calor) quer através da 
sua melhor utilização quer da sua eventual 
substituição por outros mais eficientes; 

- Introdução de tecnologias inovadoras no uso da 
energia: por ex. a BOMBA DE CALOR. 


3 Funções da Bomba de Calor num 
Edifício 

A BOMBA DE CALOR oferece a oportunidade de 
aumentar a eficiência do uso da energia, sempre que se 
pretenda dispôr de calor a temperaturas inferiores a 
100ºC, o que tem amplo âmbito de aplicação nos edifícios 
tanto no sector doméstico como nos sectores não 
domésticos. 

Com efeito, o aquecimento ambiente e a produção de 
águas quentes sanitárias e de cozinha são alguns 
exemplos de consumos que utilizam calor a 
temperaturas inferiores a 100ºC e que poderão s ser 
satisfeitos com uma BOMBA DE CALOR. 

Dado que a BOMBA DE CALOR pode ser reversível 
ela pode igualmente ser utilizada para produzir frio, isto 
é assegurar o arrefecimento ambiental. Por outro lado, 
porque a BOMBA DE CALOR é um dispositivo que 
extrai calor de uma fonte de calor (ar ou água) a mais 
baixa temperatura para outra a temperatura mais 
elevada e a sua eficiência (COP) aumenta se o desnível 
de temperatura a vencer fôr reduzido, ela está 
particularmente indicada para recuperar e reciclar 
calor. O calor rejeitado nas águas quentes sanitárias e de 
cozinha, já utilizadas, e no ar viciado rejeitado (para 
renovação de ar em edifícios com grandes taxas de 
ocupação ou exigências específicas de qualidade de ar); o 
calor libertado no condensador de grandes unidades de 
produção de frio; o calor latente resultante da 
condensação do vapor de água contido no ar extraído de 
piscinas ou cozinhas em instalações hoteleiras, 
constituem, entre outras, potenciais fontes de calor que 
poderá ser recuperado e reutilizado com o auxílio de uma 
BOMBA DE CALOR. Assim a BOMBA DE CALOR 
pode, com vantagem, desempenher as seguintes funções 
num edifício: 


- Aquecimento do ambiente e de água sanitária 
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- Condicionamento do ar (aquecimento e/ou refrigera- 
ção do ar ambiente) 

- Recuperação e reciclagem de calor degradado 

- Produção combinada de calor e frio 

- Produção combinada de calor e desumidificação 


As duas últimas aplicações da bomba de calor 
destinam-se quase exclusivamente aos sectores não 
doméstico; a opção de utilizar a bomba nos edifícios deve 
assentar num estudo dos: 


- consumos de energia pelos diversos tipos de 
utilização - valores e regimes; 

- cargas térmicas de aquecimento e de refrigeração 
para as temperaturas de utilização desejadas; 

- identificação das fontes de calor acessíveis ao 
edifício - respectivos débitos e níveis de 
temperatura. 


A viabilidade da BOMBA DE CALOR aumenta com: 


- Multiplicidade de funções que desempenhe 

- Duração de funcionamento, ou seja, taxa de 
utilização 

- Desnível de temperatura entre as duas fontes de 
calor 


4 Aplicações 


4.1 Recuperação de calor rejeitado no 
edifício para aquecimento (ambiente 
e/ou águas sanitárias) 

Suponhamos que são conhecidas as necessidades de 
energia para o aquecimento ambiente e a produção de 
águas quentes sanitárias e de cozinha. Seja Q, esse valor. 
Uma fracção X desta energia é normalmente rejeitada, 
quer através das águas quentes sanitária e de cozinha já 
utilizadas, quer através do ar viciado extraído do 
edifício. Pretende-se satisfazer as necessidades totais de 
aquecimento com uma BOMBA DE CALOR que 
recupere o calor rejeitado. 

Seja T, a temperatura (ar ou água) que se pretende 
dispôr no interior do edifício para satisfazer as 
necessidades de aquecimento; T, a temperatura a que o 
calor é rejeitado do edifício e T, a temperatura da fonte 
de calor ambiente (ar ou água) exterior ao edifício. 
Normalmente tem-se T,>T,>Ta. Seja COP, a 
eficiência de uma bomba de calor que extraisse calor do 
ambiente exterior à temperatura T, para produzir calor 
útil à temperatura T, e COP, a eficiência de uma bomba 
de calor que produzisse calor à temperatura T, 
requerida a partir do calor rejeitado do edifício à 
temperatura T,, podendo os seus valores serem 
estimados respectivamente através das seguintes 
fórmulas: 

P + 278 T + 28 
q q 


COP = 05 COP = 05 
a Tá 


(E =7)+5 (To -T)+5 
q a q r 
Note-se que COP,>COP, pois o desnível de 

temperaturas a vencer pela BOMBA DE CALOR na 

recuperação de calor é menor. À eficiência total COP, de 
uma bomba de calor que satisfaça integralmente as 
necessidades de calor e recupere o calor rejeitado 
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Q =XQ, poderá ser calculado através da seguinte 
relação 4: 


o é Ps (1) 
COP, COP, COP, 
ou seja 
T + 278 
COP =0:5 - g Cla) 


(T-T)+5-X(T -T) 
q a r a 


No caso de não haver recuperação do calor rejeitado 
(X=0), a eficiência deste sistema coincide com a 
eficiência de uma BOMBA DE CALOR que apenas 
extrai calor do ambiente exterior à temperatura T, e 
tem-se, então COP,=COP,; na situação em que se 
recupere o calor rejeitado do edifício (X >0) a eficiência 
total do sistema é sempre superior, i.e. COP,>COP,,. 

A expressão anterior está visualizada na fig.1, onde 
se representa a variação da eficiência total da BOMBA 
DE CALOR COP, com a fracção do calor recuperado X, 
para várias temperaturas do ambiente T, e do calor 
rejeitado T,, e para uma temperatura de utilização 
fixada em T,=50"C. 


08 


E 04 


X= Ur/a, 


tigura la 


Qual o efeito produzido pelo aumento da eficiência da 
BOMBA DE CALOR com recuperação de calor nos 
custos de exploração deste sistema? 

Da figura 1 ressalta o facto da eficiência da BOMBA 
DE CALOR aumentar quer com a fracção de calor 
recuperado quer com a temperatura do calor rejeitado. À 
um aumento da eficiência corresponde, para a mesma 
quantidade de calor produzida, um menor consumo de 
energia eléctrica e, consequentemente, uma redução nos 
custos de exploração. Ora o consumo de energia eléctrica 
exigido pela bomba de calor é dado por: 


4 Ver apêndice | 
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À 0.4 


0.8 
X=0 a, 
figura tb 
OQ. 
Wi = ea (2) 
á COP. 


(sendo Qear O calor útil libertado no condensador), e o 
consumo de energia primária exigido por um sistema 
térmico convencional equivalente é dado por: 


(3) 


(sendo Qrai O calor útil produzido pela caldeira e n; à 
eficiência térmica desta). Por outro lado, a poupança nos 
custos de exploração obtida com uma BOMBA DE 
CALOR, relativamente ao sistema convencional, é 
calculada do seguinte modo 5. 


COP - 6 | E sm 
a q 


T)+5 
)x 100 = ( |-26 


D(% =[ )x100 (4) 


cor, r,+ 278 
em que O representa a razão entre o custo do kWh 
eléctrico C, e o custo do kWh térmico C, 8. No caso 
presente, a poupança obtida com a BOMBA DE CALOR, 
com recuperação de calor rejeitado, relativamente a 
sistemas convencionais, resulta ser: 
(T-Ti+5- XT -T) 
a Fr a 


D,(%) =(1-20 )x100 (5) 


T + 278 
q 


Como se pode ver, a poupança conseguida neste caso 
(D,) é superior à poupança já conseguida com uma 
bomba de calor sem recuperação de calor rejeitado (D). 
Esta diferença é tanto maior quanto maior fôr a fracção 
de calor recuperado X e quanto maior fôr o desnível 


à Ver apêndice 2 

6 O custo do kWh térmico É. e dado por : (Preço de 1 kg de 
combustível) Poder calorífico do combustível x eficiência 
térmica da caldeira) 
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entre a temperatura do calor rejeitado e a temperatura 
ambiente (T,-T). 

A título de exemplo considerem-se as seguintes 
hipóteses de trabalho: 

Para satisfação das exigências totais de aquecimento 
ambiente e de produção de água quente sanitária será 
necessário dispôr de uma quantidade de calor Q, à 
temperatura de 50ºC; admita-se que cerca de 20% do 
calor rejeitado para o exterior, quer através das águas 
quentes já utilizadas quer do ar viciado extraído do 
interior do edifício, é recuperado à temperatura de 30ºC; 
considere-se que a temperatura exterior é de cerca de 
10"€. 

Para satisfazer as necessidades totais de aquecimento 
vamos considerar três soluções possiveis e confrontar os 
custos de exploração correspondentes: 


A - Bomba de calor sem recuperação de calor 
B - Bomba de calor com recuperação de calor 
C - Caldeira a gás 


No caso da BOMBA DE CALOR sem recuperação 
(X=0),o custo total de exploração vem dado por: 


= C, X Q,13.64 


e no caso da BOMBA DE CALOR com recuperação 
(X =0.20), seria: 


C,=C,xQ,/3.99 


No caso de se utilizar a caldeira a gás (por exemplo 
com eficiência 0.85) e admitindo que o custo total do 
kWh térmico C, vale, no caso do gás a granel, cerca de 
4$60, ou seja cerca de 30,6% do custo do kWh eléctrico, 
que admitimos ser 15$00, o encargo envolvido para a 
satisfação das necessidades de produção de calor vem 
dado por: 


(6a) 


(6b) 


C=C xQ/326 (7) 
É [o Ê 
Face aos valores apresentados é visível que: 


Com a BOMBA DE CALOR SEM recuperação de 
calor consegue-se uma poupança no custo de exploração 
de 10.4%, relativamente à caldeira a gás, e com 
recuperação de calor rejeitado uma poupança de 18,4%. 


Conclusão 

Os custos de exploração são reduzidos 
substancialmente com a BOMBA DE CALOR, em 
particular se esta reutilizar o calor que é rejeitado do 
edifício. 


4.2 Produção combinada de frio e calor 
A produção combinada é uma das aplicações mais 
interessantes das BOMBAS DE CALOR nos edifícios, 
em particular nos sectores residenciais não domésticos. 
Considerem-se as situações em que se pretende 
assegurar simultâneamente a refrigeração de uma zona 
de um edifício ou a produção de frio numa unidade de 
congelação e dispôr de água quente sanitária ou de 
energia térmica para o aquecimento de uma zona do 
edifício. Seja Q,ej a potência térmica exigida para 
satisfação das necessidades de frio (isto é, a potência 
térmica extraída do espaço refrigerado) e Qcai à potência 
térmica exigida para satisfação das necessidades de 
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aquecimento, ou seja, supõe-se conhecida a configuração 
da carga térmica do edifício Qrepf'Qcal- Seja Tref a 
temperatura a que se pretende a refrigeração, T, a 
temperatura a que se pretende dispôr o fluido (água ou 
ar) de aquecimento; T; a temperatura inicial do fluido 
que se pretende aquecer. 

Pretende-se satisfazer as necessidades de refrigeração 
com uma BOMBA DE CALOR e recuperar 
simultâneamente o calor libertado no seu condensador 
para assegurar, parcial ou integralmente, as 
necessidades de aquecimento. 

Supondo a BOMBA DE CALOR a funcionar como um 
sistema dimensionado para garantir as necessidades de 
frio, qual a fracção da potência calorífica (exigida pelas 
necessidades totais de aquecimento) satisfeita, em 
função da configuração da carga térmica Q,ef'Qcai do 
edifício? 

Qual o efeito do aumento da fracção da potência 
calorífica satisfeita pela BOMBA DE CALOR sobre a 
sua eficiência? 

No que se segue supor-se-á que o calor libertado no 
condensador da BOMBA DE CALOR, normalmente 
quantitativamente inferior ao calor útil necessário, será 
utilizado no pré-aquecimento do fluido; nestas condições 
a temperatura do condensador é a mais baixa possível, 0 
que optimizará a eficiência da BOMBA DE CALOR. 

A eficiência total COP, da BOMBA DE CALOR na 
produção combinada de frio e calor é dada por: 


COP, = COP (R)+COP (4) (8) 


e as eficiências de refrigeração COP. (R) e calorífica 
COP.(A) do sistema combinado são dados por: 


” =. +28 | 
COP (R) = 0.5 iz nc (9) 
(T'-=T 9+5 
q ref 
T' + 278 
COP (A)= 0.5 nas (10) 
E =T )+5 
q ref 


em que T, é a temperatura a que é libertado o calor no 
condensador da BOMBA DE CALOR - temperatura de 
pré-aquecimento - normalmente inferior à temperatura 
requerida T,. A fracção da potência calorífica satisfeita 
pelo sistema pode ser calculada do seguinte modo 7: 


e al ef SO — Lapa, ) 

“Qu Cal q dA 
dd CC o sm ser =) 
ref Qu q [ 


em que Qial é a potência calorífica libertada no 
condensador da BOMBA DE CALOR - calor de 
pré-aquecimento e Qre; é a potência de refrigeração 
exigida. 

Por outro lado, a fracção da potência fornecida é, 
também, dada por: 


Pai 
E a 
UR; 

q ' 


12 
= (12) 


7 Ver apêndice 9 
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A título de exemplo considerem-se as seguintes 
temperaturas de referência num edíficio em que se 
pretende assegurar o arrefecimento do ambiente e 
produzir água quente sanitária: 


Temperatura inicial da água de rede (T;) . 10ºC 
Temperatura da água sanitária (1) 50ºC 
Temperatura de refrigeração (T ef) TA 

e represente-se na fig. 2 a fracção da potência calorífica 


1 


08 


S 06 
Cr 
Be 
: 
S 04 
E 
0.2 
0 02 04 0.6 0.8 1 
Q 
ref/q a 
figura 2 


satisfeita pelo sistema combinado X em função da 
configuração de carga térmica Q,ef'Qcai e da eficiência 
total do sistema combinado COP.. Verifica-se que a 
partir de Qref'Qcar > 75% o sistema combinado satisfaz 
integralmente as necessidades de aquecimento (X =1) 
com uma eficiência total de cerca de 6.4. A eficiência do 
sistema combinado decresce à medida que aumenta o 
nível de satisfação das necessidades de aquecimento X, 
ou seja, à medida que a temperatura de pré-aquecimento 
T, se aproxima da temperatura de utilização desejada 
Tq: Como se vê, deverá haver um limite para a fracção 
da necessidade de aquecimento X a partir da qual a 
produção combinada deixará de ser a mais eficaz. Esta 
mesma conclusão pode ser verificada a partir do cálculo 
da dependência do consumo de energia eléctrica da 
bomba de calor com a temperatura do calor de 
pré-aquecimento libertado no seu condensador. 

Seja (We) q O consumo de energia eléctrica exigido por 
uma bomba de calor dimensionada para garantir as 
necessidades de frio à temperatura T,ef e produzir calor 
de pré-aquecimento à temperatura Te (W.), O consumo 
de energia eléctrica exigido por um tal sistema que 
libertasse calor no seu condensador à temperatura 
requerida 1, >T 

A razão entre os consumos de energia eléctrica 


UWo)'g(Welg) depende da temperatura de 
pré-aquecimento T, através da seguinte relação: 
ae ad Pag * (13) 
Wi TT + 5 
eq q ref 


e a razão entre as eficiências totais dessas bombas de 
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calor ((COPç) g(COP.);) depende da temperatura de 
pré-aquecimento T', através da seguinte relação: 


(COP, | Wj), (14) 


Pe 
A Tt y re 
(COP, (Wo, Ea +. 


E 
+ 556 

Para uma dada temperatura a que se pretende a 
refrigeração (T ,ef=-10ºC) e para uma dada temperatura 
a que se pretende dispór o fluido de aquecimento 
(T4=50ºC), a relação entre os consumos de energia 
eléctrica ou entre as eficiências totais da bomba de calor 
está representada na fig 3 em função da temperatura de 


(COP), 


figura 3 


pré-aquecimento 74, Verifica-se por exemplo que se a 
bomba de calor produzir calor de pré-aquecimento à 
temperatura de 30º€C - inferior ao nível requerido - o seu 
consumo de energia eléctrica é apenas de 70% do 
consumo de energia eléctrica que seria, se garantisse as 
mesmas necessidades de frio mas produzisse calor no 
condensador à temperatura requerida T,=50ºC, por sua 
vez a eficiência total do sistema que produzisse calor de 
pré-aquecimento à temperatura T,,=30ºC é cerca de 
1.35 a eficiência total da bomba de calor que produzisse 
calor no condensador à temperatura requerida 
T4=50ºC. 


4.2.1 Comparação da produção combinada 
de frio e calor com a produção 
independente 

Apresentou-se anteriormente a fracção da potência 
calorífica X satisfeita pelo sistema combinado em função 
da configuração da carga térmica Q,ef'Qeai do edifício 

(Eq. 11 e fig .2). Qual o valor mais interessante dessa 

fracção X* e da correspondente temperatura de 

pré-aquecimento (T',,)* por forma a garantir que haja 
benefício energético em utilizar o sistema combinado 
face à utilização de duas bombas de calor, uma 


s4 


produzindo a mesma quantidade de frio, mas a partir da 
temperatura ambiente T,, a outra assegurando o 
pré-aquecimento da água a partir da temperatura 
ambiente (em vez dele ser assegurado pelo calor de 
condensação produzido pelo sistema combinado)? 

Se COP, (R) fôr a eficiência de refrigeração do sistema 
combinado e COP (R) a eficiência de refrigeração de uma 
bomba de calor que produzisse a mesma quantidade de 
frio Qre;à mesma temperatura de refrigeração T,.f, mas 
a partir da temperatura ambiente T, (e portanto com 
maior consumo eléctrico que o sistema combinado), se 
COP (A) fôr a eficiência térmica de uma bomba de calor 
que satisfizesse o pré-aquecimento desde a temperatura 
ambiente T, até à temperatura de pré-aquecimento T, 
dados por: 


rei + TH 
COP (Rj = 05 ——— — (9) 
(o =)F 1+ 5 
qo "ref 
T + 278 
ref 
COPIR) = 0.5 (15) 
IP =P Y+E a 
ref 
TT" + 278 
CÓPIA = 05 —+ 
TM -T)+5 (16) 
q a 


então, para que o calor de condensação libertado no 
sistema combinado à temperatura de pré-aquecimento 
T, constitua um ganho energético efectivo, é necessário 
que: 


; | cora (17) 


LO = XA"O > TT Gai 
ui e Coef COPIR) COP RR). 


e a fracção interessante da potência calorífica satisfeita 
pelo sistema combinado estará limitada a : 


(17a) 


x* = Sref [—— ane ) corra) 
Qual COP (8) COPIRO: 

Se quisermos exprimir o ganho energético através da 
temperatura de pré-aquecimento mais intessante (T,)* 
que o sistema combinado deve proporcionar por forma a 
haver ganho energético efectivo, obtém-se para o seu 
limite superior, e de acordo com a equação (12): 

| (18) 
(PP =Ço=T)A*+4T 
q q l I 

A título de exemplo considerem-se os seguintes 
valores de referência de um edifício onde se pretende 
assegurar o funcionamento de uma unidade de 
congelação e produzir água quente sanitária: 


Temperatura de refrigeração (T ef) -10ºC 
Temperatura inicial da água de rede (T;)  10€ 
Temperatura da água sanitária (T,) 50ºC 


e represente-se na fig .4 
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figura 4 
(COP) (W) 
cg eq 


ai o ” tr, + 956 (14) 
(COPO, (Wo, E nd ae 


f + 556 

fracção interessante da potência calorífica satisfeita 
pelo sistema combinado (Eq. 17a) em função da 
temperatura de pré-aquecimento T., e da temperatura 
ambiente (T; = 10º€C; 30º€). 

A relação (18), ou seja, a temperatura de 
pré-aquecimento interessante (T',)* correspondente à 
fracção da potência calorífica interessante X*, está 
também representada na mesma figura. 

Note-se que às temperaturas de pré-aquecimento 
correspondem, de acordo com as relações (11 e 12), 
configurações de carga térmica igualmente 
representadas na mesma figura. 

Verifica-se, por exemplo, que para a temperatura 
ambiente T,=10ºC e para a configuração de carga 
térmica Qref'Qcal —0.58 ou para a temperatura de 
pré-aquecimento T.,=40ºC - à qual, corresponde o valor 
X=0.75 - a fracção interessante que deverá ser satisfeita 
pelo sistema combinado é X*=0.60 e o correspondente 
valor da temperatura de pré-aquecimento (T',)* =34ºC. 
Significa que, para uma dada configuração de carga 
térmica do edifício QrerfQcal, comparando os custos 
energéticos envolvidos numa bomba de calor que satisfaz 
simultâneamente as necessidades de refrigeração e de 
pré-aquecimento com os de duas bombas de calor que 
produzam o mesmo efeito (uma refrigeração e a outra em 
pré-aquecimento), é preferivel que: 


- O sistema combinado satisfaça, no máximo, 60% das 
necessidades totais de aquecimento (X*=0.60) em 
vez dos 75% que, no máximo, podiam ser satisfeitos 
face à configuração de carga térmica admitida; 
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- À temperatura de pré-aquecimento proporcionada 
pelo sistema combinado seja limitada a T,=34º€ 
em vez do valor T,=40ºC que, no máximo, podia ser 
proporcionado face à configuração de carga térmica 
admitida. 

4.3 Recuperação do calor evacuado de 
uma piscina 

Uma outra aplicação interessante da BOMBA DE 
CALOR consiste em utilizá-la para recuperar o calor 
contido no ar húmido extraído de uma piscina, o qual, 
inclui uma quantidade significativa de calor latente. 
Este poderá ser extraído se o evaporador da bomba de 
calor fôr dimensionado para produzir a condensação do 
vapor de água. Em geral, o evaporador está 
dimensionado para extrair apenas a fracção do calor 
correspondente ao calor sensível, o que, no caso vertente, 
é a parte menos significativa do calor total contido no ar 
evacuado. O calor assim recuperado com a bomba de 
calor é posteriormente reinjectado no ar interior e/ou na 
água da piscina. 

Considere-se então uma piscina semi-olímpica, em 
que se mantém a água à temperatura Tágua, O ar à 
temperatura T,; com o teor de humidade U;,, tendo uma 
utilização média diária de 12h. A taxa de renovação do 
ar exigida pela ocupação do edifício é R; a temperatura 
média do ar exterior é Tae e o seu grau de humidade 
média é U,. A taxa de calor perdida por renovação do ar é 
Q,, em que uma fracção importante deste valor (mais de 
50%), corresponde ao calor latente do vapor de água 
contido no ar extraído. 

Pretende-se utilizar a BOMBA DE CALOR que 
absorva, através do seu evaporador, o calor contido no ar 
húmido extraido da piscina Qubs (calor sensível mais o 
calor latente) de modo a arrefecê-lo até à temperatura 
exterior T,e. O calor libertado no condensador desta 
bomba de calor Qcal é posteriormente reinjectado no 
interior a uma temperatura T,. 

Supondo a bomba de calor a funcionar como um 
sistema dimensionado para arrefecer o ar húmido 
extraído da piscina até à temperatura exterior T,e € 
produzir a condensação do excesso de vapor de água nele 
contido, qual a quantidade de calor Q.ai libertado no 
condensador desta bomba de calor e qual a sua 
temperatura T? 

As eficiências de refrigeração COP (R) e calorífica 
COP(A) vêm dadas neste caso por: 


Tot 28 T + 28 
COPiR) = 0.5 COPIA) = 05 q (19) 
(T -T )+5 (T-T )I+5 
q ae q giá 
Note-se que 
= = (21) 
Real a Rio + Wes PCR Cr da) 
(COP), (Wo), T + East + 556 


e a taxa de calor libertado no condensador da bomba de 
calor. p e car São a massa volúmica e a capacidade 
térmica mássica do ar; W, é a potência eléctrica 
consumida pela BOMBA DE CALOR. 

Através das condições termohigrométricas do ar 
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interior da piscina (Tá, U,) e sabendo que o ar sai 
saturado do evaporador à temperatura do ar exterior 
(Tae, 100%), determina-se através do diagrama 
psicrométrico do ar, a variação de entalpia do ar húmido, 
ou seja, a quantidade total de calor absorvida pelo 
evaporador da bomba de calor Qabs: 

A temperatura T, do calor libertado no condensador 
dessa bomba de calor pode ser então calculada do 
seguinte modo: 


= 144 0,55 050, y, T +05 RT * pe, R a (22) 
q 0.5p Car Rh + 139p a R pio Z70 
e a quantidade de calor libertado no condensador Qual, 
que é reinjectado no interior da piscina, pode ser 
determinado através da expressão (21). 

A título de exemplo considerem-se os seguintes 
valores de referência: 


Temperatura do ar interior da piscina (T,,) 25 € 
Humidade do ar interior da piscina (U ,) 75% 
Temperatura média do ar exterior (Te) L0ºC 
Humidade média do ar exterior (U,) TO% 
Taxa de renovação do ar da piscina (R) 3000m3/h 
Taxa de calor perdido por renovação 

doar (Q,) 39kW 
Massa volúmica do ar (p) 1.15 kg/m3 


Capacidade térmica mássica (Car) 1.03 kJ/kg.ºC 

Através do diagrama psicrométrico do ar calcula-se a 
variação de entalpia do ar húmido correspondente aos 
pares de valores (25º€, 75%) e (10º€C, 100%) e que, no 
caso presente, vale cerca de 34.5 kJ/kg. À taxa de calor 
total absorvido pela bomba de calor vale então: 


Q, = 34.5 X di 
abs “3600 
e a temperatura do calor libertado no condensador é, de 
acordo com (22), T,=56ºC. 
A taxa de libertação de calor na bomba de calor vale: 


Xx1.15 =33kW 


3000 
Q =1.15xX1.03X - x (56 — 10)=45 kW 
cal 3600 


o que satisfaz integralmente a exigência de 
aquecimento, pois é superior à taxa de calor perdido por 
renovação do ar da piscina (Q =39 kW). 
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Apêndice 
L. 
E o (1) 
COP, COP, COP. 
Sabendo que: 
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A 
q, E q, * Q, 
Cy 
DB A, core (A2) 
E W+W 
a Fr 


em que 

Q,=XQ,=-W,COP, é a quantidade de calor produzido 
à temperatura T, por uma bomba de calor, a partir do 
calor rejeitado do edifício à temperatura 7 ,; 

Qu = WaCOP, é a quantidade de calor produzido à 
temperatura 7, por uma bomba de calor, « partir do 
ambiente exterior à temperatura T,. 

Da equação (Al) obtemos: 


Q,-0x, 


br as (A3) 
a COP, 
Por outro lado 
XQ 
E (A4) 
r COP, 


Substituindo as duas últimas relações (A3) e (A4) em 
(A2), obtem-se a relação pretendida. 
2. 


Deo=[ )x 100% = (1201 ) x1009(4) 
COP | 7, + 28 


Seja Ce, (COB ce Co custo de 1 kWh térmico obtido 
com um sistema convencional, o custo de 1 kWh térmico 
produzido pela bomba de calor (B.C) e o custo de 1 kWh 
eléctrico, respectivamente. À poupança relativa obtida 
com a B.C define-se como sendo: 


C- (C) 
| BC 7 
pra = (CEC ) 100% (AS) 
À | 
Mas 
C, C 
C. — D e Cao = COP 
Substituindo em (A5) vem: 
! [6] 
D(%) = |-— ) x 100% (A6) 
COP, 


Como COP=0.5 X (Ty + 278) (TT) + 5, 
substituindo em (A6), obtem-se a relação pretendida. 
3. 


pu Seat Per 88 -T, +27, ) 

“Qu A q Ay 
Fm a Pp tu8- 22 -T) 

| | Seal ' 

Sabendo que: 


Q.,.+0 A? 
COP = COP (R)+ COP (A) = Lea (AM) 
É e e W 


(A8) 


c 


Dial E V. a Jr, 
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To .+ 276 


coP (Rj=0.5 —E AS) 
| A 
q ref 
T + 278 
COP (R)=0,5 —E——— o 
i T=T 345 
q ref 
Sa (A10) 
cal F, id F, 
Combinando as equações (AT), (AB) e (A9) teremos 
o -0 (ra er LO | (A11) 
cal ÉrefA Tite 
Combinando as equações (A10) e (All) e 


rearranjando os termos, obtemos a temperatura de 
pré-aquecimento TT, em função dos outros parâmetros: 


Ê, , Q er us a 10 | 


FT as 7 Q T + 278 
T. = cal ref (A12) 
há l Dur 2 
Fr. ed cal Pa + 278 


substituindo (A12) em (Ã10) teremos finalmente o valor 
da fracção de potência calorífica X satisfeita pelo 
sistema combinado expressa em função da configuração 
da carga térmica do edifício (Q,er'Qcai) e das 
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temperaturas inicial do fluido (T;), de aquecimento 
pretendida (T,) e de refrigeração pretendida (T,ep), 
respectivamente. 


4. 
144 O bs “USA ps em + ie, RT, + ap Car R Tae(22) 
T == EU EEE 
q 0.5p Es RT e! 139 CR 
Sabendo que 
T + 278 
= ae (A 13) 
COP(R)=0.5 —E—— 
q ae 
W = O ut (A 14) 
“e COP(R) 
(A 15) 


Du WS PEÃO La 
Combinando as equações (413), (Al4) e (A15) 
teremos: 


Das (cr, — Ta) + 5) (A16) 


pe RT -T.)y= 
ar q ae 06(L io + 278) 


Cabs 


donde podemos resolver em ordem a T,, resultando, após 
rearranjo dos termos: 


“050 T tM4Q, +05pc RT +189pc RT 


T = — 
| 0.5p ba, RT + 139p Car R a 
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RESUMO 


Um problema mal formulado origina paradoxos. 
Uma vez formulado, convenientemente, estes 
paradoxos são fácilmente eliminados. Em mecânica, 
um dos conceitos geralmente considerado fundamental 
é o da reversibilidade das suas equações. Desta 
consideração deriva a bem conhecida afirmação de que 
a mecânica é irreconciliável com a termodinâmica. 
Propomo-nos mostrar que só aparentemente estes 
domínios da física são independentes e que falsos 
pressupostos estão na origem desta dificuldade de 
interpretação. 


1 Newton's law f=ma 
one of the roots of the problem 

It is common knowledge how Newton solved the 
important problem posed by the movement of comets 
and planets by introducing the well-known equation 
f=ma and the law of universal gravitation. This 
equation is prophetically valid in relativity [1] if written 
under the form f=dp/dt where p=moyv, y=(V 1-B2), 
B=vlc, my and v being, respectively, the mass and 
velocity of the body on a given frame, and c the velocity 
of light in vacuum. Problems do arise, however, when 
considering that my may vary [2]. 

The problem of the relation between mechanies and 
thermodynamics lies, in the last analysis, as will be 
seen, in the definition of quantity f. In the case of 
gravitation and of the electromagnetic field the criterion 
followed was shown to be internally consistent and was 
consecrated by experimental success. 

For gravitation or electromagnetLic forces, it is also 
possible to show, through Hamilton formalism, that 
there exists invariance in relation to symmetry in time. 

In fact, Hamilton canonical equations 


em 
and 
3H 


are symmetric in relation to time, since if t is associated 
with (-t) we have 
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ABSTRACT 


An ill-formulated problem gives rise to paradoxes. 
However, once concepts are conveniently defined, such 
false paradoxes are easily eliminated. In mechanics, 
one of the concepts generally held as fundamental is 
the reversibility of its equations. Hence, Lhe 
well-known statement that mechanics is irreconcilable 
with thermodynamies. We propose to show that only 
apparently are the two branches independent and that 
wrong assumptions are at the basis of this 
misinterpretation. 


d d d 
dt 


e ud = — 
dt  d(-t) 


Then Hamilton equations may be written 


d | je 9H 
Cage PET dq , 
d aH 
dit” dp) 
it following that 
R aH 
P, RR a Eee 
9 
and 
ai 
“agi: 
dp, 


Once we accept the fact that the domain of mechanics 
is the one of equation f= dp/dt other forces, besides those 
previously referred to, must be taken in account [3] 
which may simply result from the interaction of mass m 
with other particles [4]. This is, in fact, what happens, 
for example, in the movement of a piston confining a gas 
in a cylinder, or in the movement of a projectile in the 
air. Nothing, in such situations, is conceptually 
different, despite the difficulty of a mathematical 
definition and measurement of f. Because the piston 
moves, the value of f is altered in relation to its value if 
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it were at rest (Doppler effect), which makes it necessary 
to calculate fon the piston frame [5]. That which causes 
a false paradox to arise is, as we shall see, precisely the 
fact that the form of f, resulting from the piston 
movement, no longer allows a set of equations, like the 
ones previously referred to, which permit an invariance 
in relation to symmetry in time. 


2 The Concepts of Reversibility in 
Mechanics 

Once we introduce Lagrange function L of an 
interactive material points system, it is possible, on 
certain situations, to proceced to simplications like, for 
instance, considering two subsystems A and B and split 
Lagrangean L into Ls and Lg [6]. But nothing, in 
principle, warrants the statement according to which 
the movement of a mass m in interaction with a set of 
particles, whose number is very large, is outside the 
sphere of the most general principles of physics which 
may be called mechanical [3,7]. Now these principles 
allow to conclude that simplifications give rise to a set of 
differential equations that are not valid the whole time 
[8]. In such a case subsystem A may have a dissipative 
Hamiltonean [9] which only approximately can describe 
the interaction with subsystem B and so, both f=dp/dt 
and the equations associated the Hamiltonean cease to 
be symmetric in time, since they are incomplete |9]. 

What really happens is that a description in terms of 
dissipative forces or dissipative Hamiltoneans is not a 
complete description [10] but may be an approximate 
description for most of the time (Boltzmann, in his 
famous H theorem also introduced the hypothesis of the 
molecular chaos, obtaining a solution which is not 
complete) [11]. But in a more restricted sense equations 
with dissipative forces or dissipative Hamiltoneans 
belong in mechanics (of physics! ) in the sense that they 
are consistent with the general solution which, at least 
in principle, may be conceptually arrived at through the 
momentum-energy conservation [5]. 

Here is then a well-defined problem of reversibility. 
In fact, in a complete description, although the "system" 
[12] (for exemple, the piston in the previously described 
model) evolves towards an equilibrium point where it 
may stay for a long time (say 10º years) in which case it 
is reasonable to talk about an irreversible 
transformation (with time moving forward the piston 
and the gas will take a long time to repass the initial 
point...), the equations based on the momentum-energy 
conservation of all the entities constituing the SYSTEM 
[12] are, nevertheless, reversible. We can say, then, that 
in a given sense of reversible (sense | associated with the 
differential equation) the process is irreversible, in a 
different sense (sense Il associated with the 
momentum-energy conservation) it is reversible or 
irreversible, depending on whether one is thinking of an 
unlimited or limited time (sense Ila and IIb), and in a 
still different sense ((Ilc)-symmetry in time) it is 
reversible. 

Paradoxes will obviously and unavoidably arise if 
presuppositions are not clearly defined. 
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We shall now see that, to this difficulty, another 
problem might be added: that of the concept of 
reversibility in Thermodynamics. 

From what has been said it can be seen that the 
problem of reversibility, whether in senses IL or II, is 
related to the inicial conditions [13]. In fact, for most of 
the initial conditions, the "system" (for example, the 
piston) evolves according to the trajectory of the 
differential equation (non reversible 1), although for the 
initial conditions corresponding to the inversion of all 
the velocities of all the entities constituing the 
SYSTEM, the trajectory of the "system" is not consistent 
with the differential equation and returns to the initial 
position (reversible Id). Such conditions are therefore 
extraordinary, although mathematically possible 
[14,15]. Equations of mechanics may then be considered 
irreversible in senses | and Ilb and reversible in the 
other senses. 


3 Clausius Equation - dU=dW+dQ - 
the other side of the problem. The 
various meanings of Heat. 

Thermodynamics is generally seen as the domain of 
heat [16] defined through the so-called "First Principle 
of Thermodynamies" (FPT) dU = d&W + dQ [2]. 

But this FPT has given rise to irremovable 
contradictions [2,17,18)]. 

However, no contradictions or misinterpretations will 
arise if the Energy Conservation Principle 
(Energy-Entropy Principle) [12] is used between two 
equilibrium points by introducing quantities work (W) 
and heat (D) [12,19]: 


W=AU. 
FU=U+Up, where U is the energy of a subsystem in 
contact with a source of energy Up 


W= AU + AU ,, 
and we can then write Q=-AU pif 


AU=W+Q [2]. 
For the SYSTEM [12] 


W=AU 
and the SYSTEM is "mechanical" [20], the classification 
of“mechanical" and "thermodynamic" no longer making 
sense ( a “thermodynamic system" is a subsystem of a 
“mechanical system" !) [21]. 


3.1 The notion of equilibrium 

We can therefore say that a “mechanical System" 
tends to an equilibrium point (sense |) and, because of 
this, internal energy is a function of entropy and volume 
[22]. (Note that the definition of microscopic entropy 
associated with the maximum of thermodynamic 
probability only apparently is incompatible with the fact 
that the piston and the gas go back to the initial point in 
the phase space, since the value calculated from the 
maximum of thermodynamic probability only means 
that, most of the time, the system is to be found in the 
microstates of this maximum). 

We are thus confronted with two apparently 
contradictory equilibrium notions: a “mechanical” one, 
associated with the piston movement - the piston attains 
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equilibrium position when it "stops" at a given time (in 
fact the piston is always moving but this movement is 
physically small); the other notion is a "Lhermodynamic" 
one, associated with physical parameters which can be 
measured in the gas and possess a given mean value in à 
physically short time interval which may be replaced by 
an instantaneous quantity. This second notion is 
consistent with the first most of the time, although the 
mean value obtained for an infinitely long time 
coincides with the average obtained in short time 
intervals for most of the time intervals considered. This 
does not obviously imply that at some time intervals the 
piston does not leave its rest position ("equilibrium 
position"). In the mechanical sense previously defined 
we would say that during such time intervals we would 
have an off-equilibrium situation but, in a 
thermodynamic sense such states are equilibrium states 
- the mean values obtained for the microstates which 
maximize thermodynamic probability are the same as 
those obtained for all possible microstates. 

We may, therefore, refer to two equilibrium concepts: 

| - equilibrium in a restricted sense, related only to 

the microstates of thermodynamic maximum 
probability, generally called equilibrium microstates. 

IH - equilibrium in a broad sense, related to all 

possible microstates. 

In this second sense the non-equilibrium states are, 
obviously, the transition states between two equilibrium 
states. 

There arises, then, a notion of reversibility associated 
with dS=0 [22]. Here also, there appears in 
“thermodynamics" a terminological confusion which 
results from mixing up two concepts: the reversible 
transformation concept and the quasi-static 
transformation concept [21,23]. 

Once the incompatibility between "mechanical" and 
“thermodynamic” is proved to be false, there will be no 
incompatibility between Boltzmann and Poincaré, or 
between Carnot and Darwin [16,24]. (Note that the idea 
of irreversibility appears associated with the evolution 
of the Universe, since only as an abstraction can a 
SYSTEM be isolated). 

Perhaps the only true incompatibility is to be found in 
wrongly formulated enuntiations [22,23]; for without a 
convenient interpretation, any attempt to solve the 
“problem of irreversibility" [25] through mathematical 
simplifications that create asymmetries may only lead 
to an apparent solution for an ill-formulated problem. 


3.2 The notions of heat 
Considering what has been said, we have, between 
two equilibrium points 
W=AU . 
Ifa subsystem of energy U is assumed to be in contact 
with a source of energy Up, 


W=AU =AU+ AU 


In an infinitesimal transformation between two 
physically near equilibrium points 


dW= dU = dU + dU , 


Ifwe write 
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dU = dW + dQ 
it follows that 


dQ = — dU,. 


In sense 1, let us designate internal energy U as heat 
J 

It is obvious that, in this sense, heat of system U 
cannot coincide with U or Up (U=U+Up). This being 
so, the heat exchanged between Lhe subsystem of energy 
U and the source ofenergy Up is Q=-AU p. 

We have, therefore, a second notion of heat — heat in 
sense II, or exchanged heat which, in a module, is equal 
to the variation of heat from source AU p. 

Let us now consider the elementary variation of a 
function U (V,S) in an elementary displacement dP: 


dU = <dU,dP> 
where dU is the exterior derivative of U and <dU, dP> 
is the contraction of form dU and vector dP. 
We therefore have 


au=() av+() 
“Nav “as 


5 V 
It can be easily demonstrated that 

( 9U ) 

— | =— p 
9V /g 

and 
[ ) 
as), 


where p is the pressure and T Kelvin temperature [26]. 
Therefore 


dU=-pdV+TdsS. 
Let dU be formally divided into two (Pfaffian) 
quantities 


dW= — pdV 
and 
AQ=TdsS. 

It is obviously possible to call W and &Q respectively 
(elementary) work and heat. But to prevent any possible 
confusion with the other concepts I and II, we must 
explicitly call them heat II and work II. 

These quantities, dW and dQ, possess no energetic 
meaning that can be ascribed to them with generality (d 
is crossed precisely to indicate that we are not dealing 
with exact differentials treated with generality). 

Like previously (note that dW is not necessarily equal 


to AW) 
dW = dU +dU,, = dU —- dQ. 
In a reversible transformation dS=dS + dSp=0 and 
dW=-pdV=dU+dUp=-pdV+TdS + TdSp. 
Therefore 
dW = — pdV =dW 


aaand since 


Técnica 3/91 Técnica 3/91 Técnica 3/91 Técnica 3/91 Técnica 3/91 Técnica 3/91 Técnica 41 


Mechanics and Thermodynamics: À Semantic Problem of Classification 


dS=0=dS+dS, 


dU 
Pro dQ 
h T T 
Therefore 
dQ=TdS=M. 


In conditions of reversibility, quantity III coincides 
with quantity II. 

Also here, it is obvious that if concepts are not clearly 
defined, contradictions will inevitably arise [17,18]. 

Bearing in mind what has been said [5], let us now 
consider the previously mentioned mode] of a gas whose 
particles move perpendicularly to a piston in motion 
without friction with the walls of the cylinder 
containing the gas. 

The pressure exerted on the piston may be written 
under the form p'=p (1-2x/u) where x and u are the 
piston and particle velocity components and p is Lhe 
pressure when the piston is at rest at height x, the 
height of the piston in relation to the bottom of the 
cylinder. 

The work of this dissipative force is, in an 
infinitesimal transformation, 


AW iss = dE 
ss. 


where E represents the sum of the kinetic and potential 
energies of the piston. 

In fact, the work of the conservative force, "piston 
weight”, is 

W = — dV 
Cliria 

where V is Lhe piston potencial energy. 

The total work is, Lherefore, 


dW = EM ais + d We cai =" 
where Tis the kinetic energy of the piston, 
So, 
dW, 2“ =dE=a(m+v). 
55. 
Since 


dk + dU =0 
U being the gas energy, we can also easily write 
GU = = dk = = UM a =p dV 
where 


2x 
Pepil — -- | 


E 
as previously stated., 
[f we formally write 
U=BIV,S) 
we have 
dU=-pdV+TdS=-p'dV 
and therefore 


TdS = -(p'-p)dV. 


Now 
EV = dr:A 
where A is the piston area, 
Therefore 


p ( En já 
do “cêa 
it following that 
dS 
— = —(p'- p)jxAÃ—. 
dt gi iai 


Whether for x<0 or for x >0, dS/dt is greater than 
Zero. 

Let us consider the case where the "gas" is composed 
of one particle. Here 


2x: 
P=P == |, 
tl 


It therefore follows that 


dS DE tas À 
— = — p| — a ZA — = 
dt u F 
2%? E. 
=p im 
Since 
NKT NKT 
| = —. — O . á 
EO Ax 
dS 2% | 
Dj q 
dt ux 


assuming that for N particles the mean force is the same 
as the one calculated for one particle. 

We can thus calculate the "production" of entropy 
through an expression which is generally associated 
with reversibility: (dU =-pdV + TdS). 


4 | The Problem of Determinism 

[tis commonly affirmed that mechanies equations are 
"deterministic”", which is an entirely meaningless 
statement [27]. In fact, any presuppositions chosen 
should be clearly defined. 

There exist at least two concepts of determinism that 
must be taken into account: a scientific and a 
metaphysical one [28]. 

A certain theory associated with a system would be 
deterministic in a scientific sense if it allowed a 
knowledge of the future of that system. To this end it 
would be necessary to know the future of the Universe, 
which is obviously out of the question. Scientific 
determinism is entirely impossible, as the observer 
himself would have to be included in the Universe. 

Metaphysical determinism is linked with the idea 
that "nature knows itself” and is, in this sense, 
determined [28]. 

Now, as already seen, although in very simple 
situations, “mechanics” equation f=dp/dt allows the 
knowledge of the future (two-body problem). But 
Poincaré has shown that the introduction of a third body 
suffices to alter the idea that "“mechaniecs" is 
scientifically determined. 

In fact, as we have seen in the simplified mode! 
presented, not even f=dp/dt is exclusively mechanical, 
since it is not equivalent to the complete solution 
obtained through momentum-energy; nor is the 
complete solution scientifically determined, because 
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even for a finite number of particles the information 
which would have to be processed in order to know the 
whole future, would have to be infinite. However, the 
physical model would be metaphysically determined. 

From the previous analysis we can see that the 
problem associated with differential equations that are 
sensitive to inicial equations is not the only one that 
questions the determinism of "mechanies"; for, in fact, a 
differential equation which tends to an equilibrium 
point is not a complete solution for metaphysically 
determined physical model whose rigorous solution, 
conceivable in principle, cannot be scientifically 
determined. 


5 Conclusion 


It has been shown that it is possible to reduce 
"thermodynamics" to "mechanies", which means that 
both "mechanics" and "thermodynamies" belong to 
physics! 

The apparent impossibility of such a reduction is due 
to the fact that definitions in the domains of "mechanics" 
and "thermodynamics" should be conveniently revised: 
“mechanics" must accept dissipative forces as an 
approximation, and "thermodynamics” must accept the 
tendency to equilibrium contained in this 
approximation as part of the most general solution 
obtained through the momentum-energy conservation - 
entropy characterizes that part of energy which has a 
dissipative character. The Energy-Entropy Principle 
(E.E.P.) can replace the incorrect formulation of the 
"First" Principle of Thermodynamics (F.P.T.) and of the 
"Second" Principle of Thermodynamies (S.P.T.), since 
the E.E.P. is the Energy Conservation Principle plus the 
tendency to define a new equilibrium state. The F.P.T. is 
merely a mathematical formulation which has no 
general physical meaning, as quantities dW =-p dV and 
AQ =T dS do not allow an energetic interpretation with 
generality. The inability to understand this lies mainly 
in the fact of considering as equivalent two statements 
contained in two enuntiations that only orally are 
identical: 


[D In a “reversible" transformation, elementary 
work is dW =-p dV and heat AQ =T ds. 

H) In a reversible transformation, elementary work 
isdW=-pdV and heat dQ=T ds. 

In fact, as shown on the diagram in the appendix, in a 


transformation between two physically near 
equilibrium points, 
dW=dU = dU + dU, =dU - dQ, 
Up being the energy of a thermal reservoir. Since 
dU=-pdV+TdsS and dUp=TpdSp, 
dW=-pdV+TdS+TpdSp. In a reversible 
transformation dS + dSp=0 and dW =-pdV. Therefore 
f=Tp. This being so, dSp=-dUp/Tp=dQ/T 
Because 
dS = dS + d$,>0 
then 
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dS = — ds ,, 
da) 
AP 
dS = T 

If the transformation is reversible, dS=dQ/T and 
dW =-pdV. 

Since by definition &Q=TdS, in a reversible 
transformation dQ=%Q (only in a reversible 
transformation does &Q correspond to the heat exchange 
with the thermal reservoir). 

The use of dU=W +) as an energy conservation 
principle [17,29,30] may lead to error, as entities tW and 
dQ are confused with exchanged energy. The E .E.P. 
contains the information afforded by the Energy 
Conservation Principle and by the asymmetry between 
internal energy and work and it should not therefore be 
confused with the expression dU=&W+dQ whose 
physical significance can not be generalized. 

The various concepts of "work", "heat", "reversible”, 
“equilibrium", "mechanics" and "thermodynamics" 
originated a false problem whose solution can only be 
found upon a careful revision of the physical meaning of 
all the entities involved. 
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RESUMO 


Neste trabalho é apresentado um estudo sobre 
utilização de agregados de sensores na determinação 
de direcções de chegada de sinais. Uma revisão dos 
principais métodos é apresentada, focando com 
especial atenção os de alta resolução. À determinação 
do número de sinais incidentes também é considerada. 


| Introdução 

A recepção de sinais usando um agregado de sensores 
tem sido uma solução atractiva em problemas de 
detecção e estimação, visto que o agregado oferece a 
possibilidade de ultrapassar as limitações de 
directividade e largura de banda de um único sensor. O 
desenvolvimento dos meios computacionais velo 
permitir a utilização de técnicas de processamento de 
sinais que permitem um auto-ajuste ou adaptação do 
agregado a mudanças no conjunto de sinais recebidos. 
Nos últimos anos tem redobrado o interesse pelos 
agregados adaptativos tendo aumentado o esforço de 
investigação, de forma a conseguir reduzir a 
vulnerabilidade da recepção de sinais desejados em 
presença de interferências. São já largos os campos de 
aplicação dos agregados adaptativos que vão desde o 
Radar e Sonar à Sismologia e Geofísica, à Astronomia, à 
Medicina e, obviamente, às Telecomunicações. 
Actualmente, as principais áreas de investigação em 
agregados adaptativos são: 


a) Detecção de direcções de chegada dos sinais 
(direction ofarrival - DOA): 

b) Supressão de interferências; 

c) Obtenção de imagens com grande resolução 
angular; 

d) Conformação adaptativa de feixe e varrimento em 
situações de incerteza na localização dos sensores. 


Na figura 1 apresenta-se o esquema dum agregado, 
Como se observa, as saídas dos diferentes sensores são 
ponderadas e somadas. No caso de um agregado 
adaptativo, os pesos podem variar de acordo com as 
características dos sinais de entrada €e com o que se 
pretende que seja o sinal de saida. 

Em todos os campos de aplicação de agregados 
adaptativos, o processamento de sinais desempenha um 
papel fundamental, cabendo-lhe, entre outras missões, a 
de fornecer os meios de determinação dos pesos. Neste 
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ABSTRACT 

In this paper, a study concerning the use of arrays 
of sensors in the determination of signal directions of 
arrival. À revision of the main methods is presented 
focussing with special attention the high resolution 
ones. The determination of the number of impinging 
signals is considered too. 


figura 1 


agregado 


trabalho, vamos abordar, apenas, o caso da detecção de 


direcções de chegada. Este problema tem merecido 


grande atenção nos últimos anos, tendo motivado o 
aparecimento de várias técnicas que têm utilidade 
noutros domínios, como seja o da Análise Espectral de 
sinais. 

Supõe-se que os sensores do agregado recebem sinais 
de fontes no seu campo de visão, mas com direcções de 
chegada desconhecidas. A saída de cada sensor é uma 
soma de tais sinais e ruído. Em geral, os sinais recebidos 
pelos sensores podem ser incorrelacionados, 
parcialmente correlacionados ou completamente 
coerentes entre si. Normalmente e por razões físicas, 
supõe-se que não há correlação entre os sinais e o ruído. 
Vamos considerar sinais que são adequadamente 
caracterizados por uma única frequência - sinais de 
banda estreita - originados a grande distância do 
agregado, pelo que as ondas associadas podem ser 
consideradas ondas planas. 
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Este tema tem sido objecto de grande actividade de 
investigação nos últimos 15 anos. O impulso dado pelos 
métodos de Capon e de Máxima Entropia (1967) [5,9] foi 
de enorme influência pela motivação que introduziu na 
procura de soluções novas e mais potentes. Assim, 
partindo da análise do famoso método de Pisarenko [15] 
para a análise espectral de sinusóides em ruído, Schmidt 
[19] desenvolveu o algoritmo MUSIC (MUltiple Slgnal 
Classification) usando uma interpretação geométrica 
baseada na decomposição do espaço de observações em 2 
subespaços: o dos sinais e o do ruído. Este trabalho 
marcou profundamente a evolução futura dos métodos 
de análise espectral espacial, duma forma geral, 
inclusivamente os métodos de detecção de sinais de 
banda larga. A partir da sua publicação nunca mais 
deixou de se falar nos dois subespaços, podendo, 
inclusivé, agrupar os métodos de acordo com o subespaço 
que utilizam: 


MUSIC, MIN-NORM (norma mínima) [6] e 
PROPAGADOR [12] são formulados em termos de 
vectores do subespaço de ruído; 

- ESPRIT (Estimation of Signal Parameters Via 
Rotational Invariance Techniques) [17,18], TAM 
(Toeplitz Aproximation Method) [7] e técnicas afins 
são formuladas em termos de vectores do subespaço 
dos sinais. O algoritmo ESPRIT é extremamente 
importante, por ter uma base teórica muito simples e, 
em princípio, poder ser implementado de forma 
computacionalmente pouco custosa. O método TAM 
pode ser considerado um caso particular de ESPRIT, 
embora, a sua formulação original tenha como base a 
representação de estado dum modelo ARMA, 


No que diz respeito ao método de máxima 
verosimilhança o seu desenvolvimento tem sido 
bastante menor, dadas as dificuldades analíticas e 
computacionais que apresenta. Nos últimos tempos, 
tem-se notado um aumento de actividade na procura de 
soluções do tipo iterativo para sua implementação. 

Recentemente, foram propostas mais algumas 
alternativas interessantes que estão a ser objecto de 
investigação, nomeadamente: 


- formulação dos métodos MUSIC e MIN-NORM no 
espaço do feixe o que significa uma pré-conformação 
de feixe; 

- análise espectral espacial localizada; 

- utilização de momentos de ordem superior. 


Não faremos, aqui, a abordagem destes métodos. No 
que se refere ao estudo das propriedades estatísticas de 
todos os métodos apresentados, far-se-ão, apenas, 
referências breves ilustrando comparativamente as 
prestações dos diferentes métodos. Um estudo mais 
profundo dos diferentes métodos ultrapassa o âmbito 
deste trabalho, 

Este trabalho será constituido por três partes 
fundamentais: 


a) Estabelecimento do modelo de dados onde 
formulamos as hipóteses subjacentes aos métodos e 
algoritmos que serão desenvolvidos. 

b) Apresentação dos métodos usados na determinação 
da DOA. Consideremos os métodos Clássico, de 
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Capon, de Máxima Entropia e de Máxima 
Verosimilhança; a partir destes obteremos os 
métodos de alta resolução: MUSIC, MIN-NORM, 
ESPRIT e PROPAGADOR. Abordaremos as formas 
de os implementar, tomando como base a matriz de 
covariância. 


2 O Modelo de Dados 

Normalmente, fazem-se várias hipóteses que 
conduzem a algoritmos mais manuseáveis. Assim, O 
meio é suposto ser isotrópico e não dispersivo para que a 
propagação possa ser considerada como realizando-se em 
linha recta. Supõe-se que as fontes se encontram 
suficientemente longe para que as ondas que incidem no 
agregado possam ser consideradas planas. Suporemos 
que os sinais têm a mesma frequência central (wpy) e 
envolventes complexas que são funções lentamente 
variáveis do tempo - banda estreita. Com estas 
hipóteses, num dado ponto, cada sinal tem a forma 
s(t).eswst, onde s(t) é a envolvente complexa que se 
admite ser um processo estocástico estacionário de 
média nula e variância P, suposta definida positiva. 
Também se admite, algumas vezes, que s(t) é um sinal 
determinístico tipo potência de valor médio nulo. 
Consideraremos um agregado plano com L elementos, 
com uma geometria qualquer, onde incidem M ondas 
planas. 

Sejam 0, o ângulo entre a direcção de chegada da onda 
ie a normal ao plano da agregado!, d, a distância entre o 
sensor j e o ponto de referência? e f, (0) é a resposta 
complexa do jésimo sensor à frequência wy e direcção O, 
(ver figura 2). Tomaremos s (t).eswot como sendo a forma 
de cada sinal no sensor de referência (sensor 1), Sendo 
assim, com a ajuda da figura 2, podemos verificar que a 


frente de 


=“ &” onda 
es 

| a, 

| 
figura 2 
resposta do sensor j às M ondas é: 

M 
p= (1) 


ulg=> F(O)sid.e 
J — "| Ê l 
= 
onde 1; (0))=d,c .senQ, é o atraso entre o 1º e o jesimo 
sensores para a iºsima onda incidente no agregado vinda 


1 suporemos, apenas, à determinação do azimute e não da elevação. 
2 seguir-se-á uma prática habitual de considerar um sensor como 
referência. 
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da direcção 0,, sendo c a velocidade de propagação das 
ondas no meio considerado. 
De (1) conclui-se que a informação desejada está 
contida na representação em banda de base: 
M 
a R$ iê 
x, = 2. F, (0)s ().e 
j=1 
pelo que será, apenas com ela, que trabalharemos. 
Introduzem-se 3 vectores: 


— | 


03;68) (2) 


a) a resposta do agregado ou vector de direcção 
para a direcção 0;: 
a (0,) = Lf, (0) .e a] (3) 
fi (0) é a resposta complexa do jésimo sensor à 
frequência wy e direcção 0,. 
b) o vector das envolventes 
s(t) = Ls, (),] = À cos M] 


c) o vector do ruído, n(t), suposto independente dos 
sinais. 


(4) 


O modelo das observações vem dado por: 


x() = A(O)s(t) + n(t) (5) 
onde A (6) é uma matriz Lx M, com L>M cujas colunas 
são elementos da variedade ('manifold") do agregado. 
x(t) é um vector a tempo contínuo, pelo que deverá ser 
convenientemente amostrado. Cada um dos vectores 
x(n) chama-se "snapshot”. 

No caso particular em que 


a) os sensores são iguais e omnidireccionais ( podemos 
supor f;(0)=1) 

b) o agregado é linear e uniforme, 
1; (0) = (j-1). d. sen(0;)/c, sendo d a distância 
intersensores, os vectores de direcção assumem a 
forma 

(6) 


com tw = w Lib) 


Ú 

pelo que A(0) se torna uma matriz de 
Vandermonde. Suporemos em tudo o que se segue 
que as fontes são incorrelacionadas e incoerentes, 
pelo que o número de colunas de A (0) coincide com 
o número de fontes. 


o FP 
atbi=|le o E ae NE, du] 


O ruído supõe-se estacionário e de covariância 02 R,, 
onde R, é conhecida. Geralmente, admite-se que o ruído 
é branco ou que se efectuou uma prévia operação de 
branqueamento, pelo que R,=1. 

Os "snapshots" são acomodados nas colunas da 
chamada matriz de dados, X. Nas aplicações recursivas a 
matriz de dados é actualizada de acordo com a recursão 

EX =aX ú e 

n n-1 n 

com a e b constantes positivas menores ou igualsa le e, 
é a nésima coluna da matriz identidade nXn. Às 
constantes a e b permitem-nos manter o controle da 
influência do passado no presente e, portanto, da 
adaptabilidade do algoritmo. Nas aplicações práticas, 
não nos é possível usar um número indefinido de 
colunas. Somos, então, levados a considerar a matriz de 
dados como uma janela deslizante. Para isso, 


+ bx e 
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introduzimos uma matriz de permutação P, tal que 


| sej=i+tloul=nej=1 


Õ nos restantes casos 
Daqui resulta a formação da nova matriz de acordo 
com a recursão (ver figura 3): 


T (8) 
K 


X =axX .P+v.e 
n n-l n 


| x a 
| O — E 


esquemas de preenchimento da matriz de dados 


figura 3 


onde 
v=bx -ax, e n>L (9) 
n n | 
ce éa Késima coluna da matriz identidade. 
A matriz de covariância, R, é, na prática, estimada 
pelo produto X.XH. É fácil concluir, no caso de uma 
janela de dimensão indefinida, que: 
R =aR + Bx x” a 
n n=] n n 
coma=a2e B=b2. No caso de uma janela finita: 


H H 
E ax,.xX, 


A possibilidade de controlar a influência do passado 
no presente através das constantes a e 5, faz com que o 
esquema da eq. (11) não tenha nenhuma vantagem sobre 
(10), pelo que a alternativa se põe entre (8) e (10). Em 
termos estatísticos o esquema da eg. (10) é melhor que o 
traduzido em (8), devido à efectuação de uma operação 
de média. Nas aplicações recursivas o uso da covariância 
conduz a algoritmos melhores do ponto de vista 
numérico [18]. Em tudo o que se segue, sempre que 
falemos em matriz de dados, supô-la-emos gerada de 
acordo com a eq. (8). De forma análoga, suporemos que a 
autocovariância é gerada de acordo com (10)3, 

Todas as técnicas que seguidamente se apresentam 
têm aplicação imediata na análise de séries temporais, 
desde que consideradas como sendo a soma de sinusóides 
moduladas em amplitude por sinais de banda estreita. 
Na figura 4 mostra-se a forma de gerar um dado 
"snapshot" para tal situação. 


o (11) 
R. = GR. + Dx x 


x (n) = [ x () x (m) x) x (m) em ] 


figura4 “snapshot” de serie temporal 


3 Métodos de Determinação de DOA 


ul Técnicas convencionais 


Como se disse atrás, um agregado pode ser 
considerado como um sistema com L entradas e uma 


3 para acovariância cruzada a expressão é análoga. 
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saida. Esta obtem-se pesando cada sinal de entrada e 
somando [12,14]: 
o) H (12) 
Y(in)=w .xtn) 

onde w é o vector dos pesos, funções de O e x(n) é o 
"snapshot", A escolha de w depende dos sinais de 
entrada e do sinal desejado; normalmente, escolhe-se w 
de acordo com um dado critério que tem como base a 
extremização dum dado funcional. Por exemplo, w pode 
ser escolhido de forma que a potência de saída seja 
maximizada. Um outro critério muito útil consiste em 
maximizar a relação sinal-ruído. 

Vamos apresentar, em seguida, os métodos mais 
importantes de determinação de DOA. 


3.1.1. O conformador de feixe 
Suponhamos que pretendemos "focar" uma dada 
direcção 6;, que, sem perda de generalidade, vamos 
supor coincidir com a DOA da onda 1, Para isso, 
escolhemos w igual ao vector de direcção a(0,). 
Consideremos que os sensores são omnidireccionais e 
iguais, com f(0,)= 1, e que não há ruído. Obtemos, assim, 
à saída do agregado um sinal y(t): 
M a 
y (8) =L.s (0) + > s (8. 3. e 
1=2 j=1 
que mostra, que a menos que a sobreposição dos sinais 
não desejados dê um sinal nulo, não podemos obter o 
sinal desejado livre de interferências. No que se refere à 
potência, os cálculos são fastidiosos, pelo que não os 
faremos aqui. Pode mostrar-se que há um máximo da 
potência na direcção da fonte. 
No caso geral de uma direcção qualquer e nas 
condições gerais indicadas na sec. 2, a potência de saida 
pode escrever-se: 


— =t 48.) 
ice da 0) k: 1 


] 
O, = a/(0) R. a(0) do 
Tomam-se os valores de O para os quais Pçcp tem 
máximos como sendo as DOA, À este estimador dá-se o 
nome de conformador de feixe [14]. No caso de uma 
única fonte, P mede a potência efectiva quando 0 
coincide com a direcção de chegada. No caso de várias 
fontes isto já não acontece, mesmo se as fontes forem 
incorrelacionadas; com efeito, os picos afastam-se das 
posições correctas e, se estiverem “demasiado perto”, 
podem fundir-se. Isto mesmo na ausência de ruído; em 
ruído,a situação piora, pelo que este estimador não tem 
utilidade. 


3.1.2 O Método de Capon 

Os problemas levantados pelo estimador anterior 
levaram Capon [3] a reformular o problema nos 
seguintes termos: 

Calcule-se w de forma a que, para uma dada direcção, 
sejam minimizadas as contribuições vindas de outras 
direcções: 

min (14) 
Lu 

Este é chamado estimador de variância mínima. 
Trata-se de uma minimização com restrições, que 
efectuada conduz a: 


w! R.w sujeitoa w! a(0)=1 


ai 
R a (15) 
ne H o-l 
ak a 
Desta forma, a potência de saída é dada por: 
l 
PT TRT, si 
a R a 


Um raciocínio muito simples permite-nos mostrar que 
dB O 


-10 


= 20 


figuraS  Conformador de feixe vs metodo de Capon. 


este estimador tem uma maior resolução do que o 
conformador de feixe. Com efeito, seja 


l 
| = — a! Rº Ra 
E, 


Usando a desigualdade de Schwartz conclui-se, de 
imediato, que aiRa.a!R-lazL e, portanto: 
Pcp=L Pym. Esta maior resolução é responsável pelo 
êxito inegável que este estimador teve. Na figura 5 
ilustra-se este facto. 


3.1.3 O Método de Máxima Entropia 

O éxito conseguido pelo Método de Máxima Entropia 
(MEM) em Análise Espectral levou a que ele fosse 
aplicado na determinação de DOA's. Neste método, 
[5,9,14] minimiza-se a potência do erro associado à 
predição da saida dum dado sensor à custa das restantes 
L-1 saidas. O problema pode formular-se, também, em 
termos de uma minimização da potência de saída com 
restrição. Escolhendo o primeiro sensor como referência, 
temos: 

Calcule-se w de forma a que seja minima a potência do 
erro de predição da saída do 1º sensor em termos das 
restantes saidas 

e w” R.w sujeitoa wºe = | ea 
ondeepéalê2 coluna da matriz identidade. 

Procedendo à minimização de forma análoga à 
anterior, obtem-se o vector w cujos coeficientes são os 
parâmetros do filtro preditor: 


l 
To-l 
e, R e, 
A “resposta angular” desse filtro permite-nos definir 
um estimador da densidade angular de potência: 


Ww = Re (18) 


| 


Tam = 1 
e,R e, (19) 
ja” R'!e [É 


Dye = 
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Este estimador tem maior resolução que os 
anteriores. Com efeito, o estimador VM pode ser 
considerado como uma “associação em paralelo" (média 
harmónica) de L estimadores ME de ordens 
sucessivamente crescentes, donde a sua menor 
resolução. Para deduzir este facto basta relembrar que a 
inversa de R tem uma decomposição do tipo Cholesky 
R1=T.TH. A matriz T é uma triangular inferior sendo 
as suas colunas os coeficientes do filtro preditor de 
ordens sucessivamente decrescentes, pelo que podemos 
obter, quase imediatamente: 


l ] ] 
a cat Eloa Ao A. 


D,m La 


sendo D'yg o estimador ME de ordem i (correspondente 
ao uso dos primeiros i sensores). 

Porém, é conveniente não esquecer que o estimador de 
Capon mede a potência enquanto que o de ME mede a 
densidade de potência pelo que a expressão anterior 
sugere um efeito de integração. Este efeito é o 
responsável pela menor resolução, mas também pela 
grande estabilidade estatística do estimador VM. Assim, 
verifica-se que o estimador ME tem a maior resolução 
dos 3 métodos apresentados, enquanto que o estimador 
VM tem a menor variância (ver figura seguinte). 


! 
Dye 


-90 -60 -30 0 30 60 90 
figura 6 Método de Capon vs MEM. 


Estes métodos têm todos em comum a particularidade 
de suporem um modelo aproximado (MA para o CF e AR 
para os outros) quando o modelo correcto é o modelo 
ARMA, 

Do estimador de Capon, pode obter-se um estimador 
de densidade de potência. Para isso, pode-se dividir Pym 
pela norma quadrática do vector w: 


” Dar a! RT | a 
a” Ra 


BL Ee ag 
A este estimador dá-se o nome de estimador de 
Borgioti-Lagunas [11]. Na figura imediatamente a 
seguir ilustra-se o comportamento deste estimador. 


(20) 


3.1.4 O Método de Máxima Verosimilhança 

Como em muitos outros problemas de estimação de 
parâmetros, o método de Máxima Verosimilhança 
(Maximum Likelihood - ML) [10] fornece-nos o 
estimador óptimo. Contudo, as dificuldades inerentes a 
este tipo de métodos são tais que o tornam difícil de 
implementar e, muitas vezes, ineficiente, obrigando a 
aproximações que o tornam sub-óptimo. No caso 
presente, embora se conheçam algoritmos que conduzem 
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= -b0 =J0 0 30 so ao 


figura 7 Ilustração da prestação do estimador de Borgioti-Lagunas. 


à estimação óptima, continua a ser bastante dificil o uso 
deste método, dado o seu peso computacional. 

Consoante se adopte um modelo determinístico ou 
estocástico para os sinais, assim se obtem um método ML 
diferente. Por exemplo, no caso determinístico, a função 
densidade de probabilidade para as observações 
escreve-se: 


— NL N 


fiX,8,8S)=n (det(RAD x 
N (21) 


q H ol 
=. [xtk) — ALOstki] RA, Jxtk) — ACO) stko] 


l 
Xe 


onde X, O e S, são respectivamente, a matriz de dados, o 
vector de direcções e o do sinal e N o número de 
“snapshots”". A maximização da expressão em (20) é 
equivalente ao problema de mínimos quadrados não 
lineares: 

Para uma dada matriz X, minimize-se 


N 
(0, S;X) = > Ix(k)— A(O). stk)P! Lx) — A(O) (he) |122) 


| 
em relaçãoa des. 

Os ângulos e sinais óptimos formam a melhor 
aproximação aos dados, no sentido em que o erro 
quadrático é minimizado. A dupla estimação de sinais e 
ângulos pode ser desacoplada e efectuada em 2 passos: 

1 estimação de sinais 
Rr (23) 

s(k)=A x(k) 

onde A-l é a pseudo inversa de A. 
2 substituição de (23) em (22), de forma a que a 

minimização em 0 se reduza, depois de se fazer 

R=I/N.X.XH,a: 


O=FMer (A!R.AI 


ou seja, o óptimo O maximiza a norma da projecção 
dos dados sobre o chamado espaço dos sinais. Não é 
difícil concluir que (23) se pode escrever na forma: 


(24) 


M 
O = ii >. la; R. a,) (25) 


1=1 
sendo a; e ay! as linhas de ordem i de A e da sua 
pseudo-inversa. 
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A solução de (24) ou (25) deve ser procurada 
iterativamente, para o que foram apresentados vários 
algoritmos sendo de referir o algoritmo EM proposto por 
Miller e Fuhrmann [10]. Não nos debruçaremos mais 
sobre este método, que dado o seu peso computacional 
não é útil em situações de movimentos de fontes. Por 
outro lado, os algoritmos subóptimos que 
apresentaremos de seguida são muito simples e 
fornecem resultados satisfatórios. De referir que no caso 
de uma única fonte o conformador de feixe é um 
estimador ML. 


3.2 Técnicas baseadas na decomposição 


em subespaços 


3.2.1 A decomposição em subespaços 
O modelo de dados antes apresentado conduz ao 
seguinte valor expectável da matriz de autocovariância: 
NE H 2 (26) 
R = kIR] = A(O)R A (0) + 0“1 
onde R, é valor expectável da autocovariância de s(t), 
suposta definida positiva Na hipótese de sinais 
incorrelacionados, é uma matriz diagonal cujos 
elementos não-nulos são as potências dos diferentes 
sinais. O objectivo último do processamento em que 
estamos interessados é encontrar A(0), Como objectivo 
intermédio pretende-se encontrar M vectores 
linearmente independentes contidos no subespaço 
gerado pelos vectores de direcção: subespaço de sinal. 
No caso de ausência de ruído, aqueles vectores são os M 
autovectores correspondentes aos autovalores não nulos 
da matriz de autocovariância, os restantes geram um 
subespaço ortogonal ao subespaço de sinais, No caso 
ruidoso, os M autovectores correspondentes aos M 
malores autovalores geram o subespaço de sinal, Os 
restantes autovectores geram o subespaço de ruído. 
Portanto, a matriz R tem a seguinte autodecomposição 
(AD): 


| (27) 
R=VAV! 
sendo os valores próprios dados por: 
n, + o k=1,..,M 
A, = (28) 
o? k=M+1,...,L 


onde os up são os autovalores da matriz A(0) RGAH(0) e a 
matriz V pode ser particionada em V=[V.: Val, tendo 
as duas submatrizes significados evidentes. 


3.2.2 Matriz de dados vs matriz de 
covariância 

Os subespaços podem ser calculados quer pelo recurso 
à autodecomposição de KR, quer pela decomposição em 
vectores singulares (SVD) da matriz de dados de X. Os M 
vectores singulares correspondentes aos M maiores 
valores singulares formam uma matriz de característica 
M, que melhor aproxima X no sentido de mínimos 
quadrados, 

É fácil ver que a AD de R e SVD de X dão a mesma 
estimativa dos subespaços, de facto, pondo X= VLZUIH, 
R=X.XH= VE2VH, 


A maior dificuldade no uso da AD de R reside na 
multiplicação matricial, dado que os algoritmos usuais 
necessitam que se forme R explicitamente. Contudo, ao 
formar R estamos a efecluar uma média o que Lraz 
vantagens do ponto de vista estatístico. À comparação de 
(8) e (10) pode levar-nos à conclusão de que depende da 
aplicação prática particular a opção de um dos 
esquemas. Num primeiro relance, (8) parece ser melhor 
em cenários de variação relativamente rápida e (10) nos 
restantes casos. Contudo, o equilíbrio que podemos 
estabelecer entre o passado e o presente através das 
constantes a e À torna os dois algoritmos equivalentes 
deste ponto de vista. De qualquer modo frise-se, que, O 
uso de (10) será preferido se se puder evitar a 
multiplicação matricial, ou se, como sucede em 
Astronomia ou Espectroscopia, o sinal não for 
directamente acessível, mas sim a sua covariância. Os 
métodos apresentados nas secções seguintes conduzem a 
valores correctos quando se usa a autocovariância 
exacta, independentemente da relação sinal/ruído e da 
separação angular das fontes. Contudo, quando se 
estimam as autocovariâncias a partir dum número finito 
de "snapshots”, aquelas técnicas exibem descentragens 
dependentes da relação sinal/ruído, da separação 
angular das fontes e, até, das especificações do agregado. 
No entanto, é possível aumentar as prestações do 
estimador no caso de um agregado linear uniforme. Tal 
melhoria consegue-se à custa da introdução do 
“snapshot” reverso ("backward") [81]. 


3.2.3. MUSIC e MIN-NORM 

A noção de subespaço vai permitir-nos a introdução 
do algoritmo MUSIC, que resulta, directamente, da 
decomposição que acabámos de efectuar na matriz de 
covariância. Essa decomposição permite que possamos 
escrever: 


RV =0".Y 


mas, atendendo a (20), temos: 


n 


A(O)R, AO) +01 V, =0".V 


n 
donde se tira fácilmente 
| (29) 
A“(0).V =0 
ou seja: Os autovectores associados ao subespaço de ruído 
são ortogonais aos vectores de direcção correspondentes 
às DUA 's efectivas. 
Assim pois, para um dado a(0) e sendo v;(i=1,...) os 
vectores próprios, os picos da função: 


MUSIC 
MUSIC 6, 


2 
> ta!y (30) 


p 


di 


1=M+4] 


ou, equivalentemente, os zeros de 


LR e 
> ja!!y j 

1=M+1 

corresponderão às direcções correctas. De notar que (30) 
não mede a potência, como nos casos anteriores, 
servindo, apenas, como detector. Este algoritmo 
apresentado por Schmidt [18] em 1979 constitui uma 


(31) 
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generalização do algoritmo de Pisarenko [15] usado em 
análise espectral de sinusóides em ruido e que usa, 
apenas, um autovector do subespaço de ruído. 

Uma alternativa à dedução atrás feita consiste em 
aproximar R-1 pela matriz que se obtém por eliminação 
dos autovectores correspondentes ao subespaço do sinal. 
Esta situação corresponde a supor que os sinais têm uma 
potência infinita, o que explica a superresolução deste 
método. Procedendo a tal aproximação na expressão que 
define o estimador VM, este transforma-se num detector 
angular que pode escrever-se na forma: 


NUS apa —— 
Lo. 9 | 
s 1! ay | (32) 
=M+1 
Se admitirmos que A,=02 (i=M+1,..., L), obtemos o 


estimador MUSIC, de que (32) é a versão ponderada 
proposta por Johnson [4]. Esta versão comporta-se 
melhor que a normal, sobretudo, em situações de relação 
sinal-ruído baixa. 

Duma forma inteiramente análoga a esta podemos 
obter o método MIN-NORM (norma mínima) pesado. 
Basta fazer a aproximação de R-1 em (18), obtendo-se, 
imediatamente: 


l 
pn DEM (a) q 


L 
| ) Nav ve 
| | EE 
r=M+1 


* qa 


-90 -60 -30 à) 30 60 90 


figura 8 MUSIC comparado com ME e VM 
Fazendo, como atrás ÀAj=02 (i=M+1, ..., L), obtemos o 
estimador MIN-NORM: 
MIN-NORM l 
pb ta)= 
= (34) 
| > a WYY oi 
1=M+1 


Este estimador foi proposto por Kumaresan e Tufts [6] 
e introduzido com base em argumentos ligados à 
predição linear. À ideia básica é a seguinte: 

Construa-se um vector d contido no subespaço de 
ruído e tal que o seu primeiro elemento seja unitário e 
tenha norma mínima. No caso dum agregado uniforme, o 
polinómio cujos coeficientes são as componentes de d, 
tem os zeros correspondentes às DOA's sobre a 
circunferência unitária e os restantes no interior dela e 
com uma distribuição aproximadamente uniforme. Pode 
mostrar-se [6] que, à parte um factor, 

(35.a) 


d=V.u 
n 
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onde u é o conjugado do transposto da 1º linha de V, pelo 
que (34) se pode escrever na forma: 


MIN-NORM l 
pb ta)= 


3 (35.b) 
[aa 

Como seria de esperar, o método MIN-NOURM está 
para MUSIC, assim como o estimador ME está para VM. 


Com efeito, podemos afirmar que [14]: 


a) O estimador MIN-NORM tem melhor resolução do 
que MUSIC; 

b) MUSIC é mais estável do ponto de vista estatístico 
do que MIN-NORM; 

c) Existe, também, uma relação do tipo "associação 
em paralelo", idêntica à descrita em 3.1.3, 
relacionando MUSIC e MIN-NORM. 


3.2.4 O Propagador 

Os algoritmos MUSIC e MIN-NORM exploram, de 
uma forma inteligente a estrutura da matriz de 
covariância, com a introdução e aproveitamento da 
noção de subespaço, gerado por um conjunto de vectores 
próprios. Contudo, cada subespaço pode ser gerado por 
outra base qualquer. No caso do subespaço de sinal, o 
conjunto dos vectores de direcção são uma alternativa 
válida. O Propagador resulta do facto de a matriz de 
vectores de direcção ser uma matriz determinística de 
característica M [12]. Nesta situação, é possível 
exprimir as restantes L-M linhas em termos das M 
primeiras 4. Sejam 


'y 


0 
Im 
duas matrizes de dimensões MXL, sendo Iy e Im 
matrizes identidades. Sejam, agora, A, e A; duas 


submatrizes de A dadas por: 


Ee (36) 


A e J 


Pei 


” (37) 


est" E 
A =) A e A=J, gy A 


Comecemos por notar que, em geral, A, é invertível, 
pelo que existe uma matriz W, (L-M) x M, tal que 
(38) 
WE.A =A, 

A matriz W chama-se Propagador. Reescrevamos (38) 
na forma: 
(39) 


=1Wº.-IA=0 


A 
H.. 5 
[W CT] A. 


L 
Esta equação é muito importante, porque mostra que os 
vectores colunas de 

W 

E) 


são ortogonais aos vectores de direcção, pelo que 
pertencem ao subespaço de ruido, podendo, então, 
substituir os vectores próprios nos algoritmos MUSIC ou 
MIN-NORM. 

Suponhamos que não existe ruído e que retomamos a 
eq. (5). Obtemos fácilmente: 


* Eventualmente, podemos considerar outras submatrizes. No caso 
especial de um agregado linear uniforme, é útil, como se verá 
mais tarde, tomar as duas submatrizes Lx(L-1) que se obtêm 
suprimindo a 1º e última linhas de A. 
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[W. -NAS=[Wº.-I)X=0 o) 
onde S é uma matriz cujas colunas são sucessivos 
vectores de envolventes e X a correspondente matriz de 
dados suposta particionada na forma 


X. 
5 


X 
Ê d 


X= (41) 


É fácil obter, agora: 

wº. X =X, 
Por outro lado, sendo a autocovariância, R=XXH, 
usando (41) e (42) obtemos: 


M x! X x! wi! 
5 5 8 


WX x! WX ww! 
E Ss Ss 5 
que nos demonstra que W pode, alternativamente, ser 


calculada a partir de. 


W R. o R (44.a) 


(42) 


(43) 


RWº-=R (44.b) 
e d 
onde ie s têm os significados anteriores e d e e significam 
direito e esquerdo, respectivamente. 
No caso de haver ruido, a igualdade (42) não se 
verifica, havendo uma matriz de erro: 


E=x-Wwºx 
Ê É] 


Esta equação mostra que W é um preditor multicanal, 
generalização de (18). 

A equação (45) é a base da estimação de W. A 
minimização da potência de cada um dos erros definidos 
por (45), ou seja, do traço de E. EH conduz a: 


(45) 


R=WR o 
Ls ss 
onde R;, é a covariância cruzada entre as partes inferior 
e superior do snapshot e Rs, é a autocovariância da parte 
superior, sendo estimadas, de acordo com o que se disse 
anteriormente por: 
R.=X.X] «e Bo=x.x 

LS ! & su b& % 

Outra possibilidade está em trabalhar directamente 
com a matriz de dados X. Neste caso (45) pode 
escrever-se na forma: 

x, = Wºx a 


A solução de norma minima é: 


É 


Wi = X | x (48) 
onde X, ' é a pseudo-inversa de X.. 

(Qualquer das relações (44.a e b), (46) ou (47) pode ser 
usada na determinação de W. Para obter resultados 
mais fiáveis, é conveniente efectuar a diagonalização 
simultânea das duas matrizes envolvidas e tomar as 
componentes principais de cada uma delas. 


3.2.4 Técnicas de rotação de Subespaço 
3.2.4.1 ESPRIT 


Há várias formulações correspondentes às técnicas de 
rotação de subespaço. A mais conhecida de todas é a 
ESPRIT [17, 18] que pode ser implementada de várias 


dB 
+ 10 
Ô 
- 1Q nr | 
- 90 = GU =50 O 30 bQ 90 
figura 9 Ilustração da prestação do propagador 


formas. As formulações mais vulgares dos algoritmos 
ESPRIT baseiam-se na consideração de que o agregado 
possui uma invariância de deslocamento, ou que é 
constituido por L/2 dobletes (ver fig. 10). 

Seja P uma matriz diagonal cujos elementos são 
elementos da forma eJ2n(dicisenb onde O é a direcção de 
chegada (DOA) e d a separação entre sensores de um 
mesmo doblete (deslocamento). Então 


d | 
— jên=sen 0. 


P = diag|e 


(49) 
E=1, co 
Sejam x e y os "snapshots”" obtidos a partir de cada 
conjunto de sensores e Ay e Ay as correspondentes 
matrizes de direcção que se relacionam através da 
matriz 

Y x 
Por outro lado, atendendo a (38) vem: 
WA-AM po 

I x 

donde se conclui que os vectores colunas de A, são os 
vectores próprios do propagador correspondentes aos 
valores próprios eJ2nd/c)senb,, 


1 2 a 4 5 [5 7 
-— o — = — e pe eo O — —— — — — 
No 
N 1 ê 3 4 5 6 7 
ARNO bn — .— — — 


figura 10 Agregado “tipo ESPIRIT” 


No caso particular de 2 subagregados uniformes, A, é 
uma matriz de Vandermonde e W é uma matriz na forma 
companheira. 

Se atendermos (50), os "snapshots" x e y escrevem-se 
na forma 


I=As e = AD. (52.2) 


donde, de forma análoga, as correspondentes matrizes de 
dados são: 
L=A8 e Y=AO,S 
Es X 


Isto significa que, atendendo a (51), X e Y se 
relacionam por: 


(52.b) 


(53) 
Ye w”x 
de forma idêntica à encontrada no propagador. 
Introduzamos as matrizes de covariância R e Q 
(54) 


R=XX” e q=xy” 
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Os respectivos valores espectáveis R e Q vêm dados, 
respectivamente, por: 
(55) 
R=A PA! 

ps A 

(56) 
Q=A PQ AÍ 

Consequentemente, valores próprios do "pencil" RAQ 
são os elementos diagonais de O: 


. . (57) 
R-AQ=APA!-NAPO AÉ=AP(I-AOD A! 


Concluimos que, para determinar os DOA, basta 
calcular os valores próprios associados ao "pencil" RAQ. 
Temos, então, 2 hipóteses: 


| Cálculo de W e determinação dos seus valores 
próprios 
2 Determinação dos valores próprios do "pencil" (57) 


O número de algoritmos para a determinação de DOA 
a partir de (53) é enorme. Não nos debruçaremos sobre 
eles pelas razões antes apontadas. 

Consideremos, agora, a situação em que o sinal tem 
ruído. Em [17] propõem-se a solução dada pelos valores 
próprios associados ao “pencil" (R-021)-1Q, sendo 02 a 
potência do ruído. Os autovectores da matriz da 
autocovariância R permitem, teoricamente, tornar a 
matriz Q triangular superior, o que não acontece na 
prática. Tal facto conduz a maus resultados. Para 
resolver o problema foram propostas duas alternativas 
que têm em comum uma manipulação simultânea das 
duas matrizes. O algoritmo mais simples foi proposto por 
Roy e Kailath [18] e consiste no uso dum critério baseado 
nos mínimos quadrados totais no cálculo da matriz !F. 
Uma alternativa a este esquema chama-se GEESE 
(GEneralized Eigenvalues utilizing Signal Subspace 
Eigenvectors) [14] e utiliza directamente os vectores 
próprios do subespaço de sinal. Outras alternativas 
baseiam-se em (53). Como já referimos atrás, o uso de 
métodos baseados na covariância conduzem a melhores 
resultados e dão a possibilidade de implementar 
melhores algoritmos recursivos, pelo que não nos 
debruçaremos sobre os métodos baseados na matriz de 
dados. 


3.2.4.2 ESPRITeTAM 


Merece atenção especial o caso em que em vez de um 
agregado de dobletes, temos um agregado linear 
uniforme. Uma forma de aplicar o esquema ESPRIT a 
este caso consiste em tomar dois subagregados com L-K 
sensores (lalL-Kek+lal) 


subagregado 1 
qui — —  — — — e 
1 2 E 4 5 É 7 
---—r.—-—-—+ — e. — — e — - 
| ágei — subagregado 2 Fes + 


figura 11 subagregados de um agregado linear uniforme 


O caso K=L/2 coincide com o discutido na secção 
anterior. Se K <L/2, então é fácil mostrar que há uma 
pequena diferença relativamente ao discutido antes, na 
medida em que o valor médio da covariância cruzada 
vale agora 
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E o (58) 

Q = A POD A. + 0 J 
sendo JK uma matriz cujos únicos valores não nulos 
valem 1 e estão sobre a subdiagonal secundária K. 

Se O for a matriz definida na secção anterior, os 
valores próprios do "pencil" 

(R-021)-M(Q-0"J,.) em 
são estimativas de DK, 

O caso particular em que K=1 é equivalente ao 
algoritmo TAM (Toeplitz Aproximation Method) [7] que 
se obtém a partir do modelo de estado para sinusóides em 
ruído. 

Sejam F a matriz da dinâmica e z(n) o estado no 
instante n, então 


nlnd = Pati (60) 


ct z(0) = Zo 


sendo F a matriz da dinâmica cujos valores próprios são 
eJtnidicisenbdie d agora a distância entre sensores. À saída 
do sensor j º é 


x, (mn) = H.z(n—-j+1) (61) 
donde se deduz que: 
f (62) 
atm= RAF. Fa=H.6.4, 
e a autocovariância 
(63) 


R=H.G.2,.6".H” 


com Zo= Elzo.zo |. A eq. (63) mostra que a matriz de 
covariância é factorizável na forma 


R=-6. 1” 
e, atendendo à forma particular de cada factor, podemos 
obter a matriz F a partir de qualquer deles. Assim, se O; 
e 95 forem duas matrizes obtidas de O por supressão da 
1ºe última linhas, respectivamente 

Fo, á 9, 

Em termos da matriz de direcções, F coincidiria com a 
matriz P e estariamos na situação em que 
consideraríamos dois subagregados com L-1 sensores, 
pelo que o método TAM pode ser considerado um caso 
particular de ESPRIT e, atendendo à forma de cálculo, 
idêntico ao algoritmo GEESE. 

No que se refere às prestações dos algoritmos do tipo 
ESPRIT, o resultado mais interessante especifica que o 
quociente entre as variâncias assintóticas dos DOA's 
detectados por ESPRIT e MUSIC é superior a 1 [20], 
donde se conclui que o algoritmo MUSIC é melhor. 


(64) 


(65) 


3.2.6 Determinação do número de fontes 
Até agora, procedeu-se como se o número de fontes 
fosse conhecido “a priori”. Na prática, tal número, M, é 
em geral desconhecido, tendo de estimá-lo a partir dos 
dados. No entanto, nas condições apresentadas na sec. 2, 
podemos explorar a estrutura da matriz de covariância 
para estimar M. Assim, se tivessemos acesso ao valor 
expectável da covariância, o problema era trivial, visto 
aquela ter L-M valores próprios iguais. Contudo, tendo 
apenas uma estimativa da covariância, o problema deixa 
de ser trivial, visto que, com probabilidade 1, os 


5 Toma-se como referência o sensor |. 
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autovalores serão todos diferentes. Para solucionar a 
questão, foram propostos testes baseados na 
característica assintótica gaussiana dos valores 
próprios. Estes testes usam o critério da razão de 
verosimilhança para a igualdade dos L-M menores 
valores próprios e que é igual à razão entre as médias 
aritmética e geométrica daqueles valores próprios [14]: 


(66) 


A partir desta grandeza obtêm-se testes que nos dão, 
de imediato o valor de M. Os mais usados são o AIC 
(Akaike's Information Criterion) e o MDL (Minimum 
Description Length) [14]. Ambos seleccionam o valor de 
M que conduz a um mínimo de uma função. Essa função 
define-se para AIC e MDL, respectivamente, por: 


(67) 
+4LM-2Mº +92 


AIC(M)= -2N.log | r y(A) 


( 
MDL(M) = —2N(L-M) . logl r MI] + [2M -M2+1].10g N 


sendo N o número de “snapshots" usados. 
Estes dois critérios têm sido usados sistematicamente 
com resultados satisfatórios. 


5 Conclusões 

Neste trabalho fez-se uma revisão geral dos mais 
importantes métodos de determinação de direcções de 
chegada. Por questões de espaço não se abordaram os 
métodos de estimação localizada e no espaço do feixe 
[1,2] que são objecto de investigação muito activa, assim 
como dos métodos baseados em momentos de ordem 
superior. À partir do que se mostrou, podemos concluir 
que os métodos de alta resolução, embora mais pesados 
computacionalmente, estão baseados numa formulação 
relativamente simples e que, paralelamente, a boas 
prestações têm uma fácil implementação. Estudámos 
essa implementação e mostrámos que os algoritmos com 
base na covariância são mais estáveis do ponto de vista 
numérico/estatístico. Finalmente abordámos a questão 
da determinação do número de sinais para o que 
apresentámos dois critérios. 
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